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We consider serious conceptual problems with the application of standard perturbation theory, 
in its zero temperature version, to the computation of the dressed Fermi surface for an interacting 
electronic system. In order to overcome these difficulties, we set up a variational approach which is 
shown to be equivalent to the renormalized perturbation theory where the dressed Fermi surface is 
fixed by recursively computed counterterms. The physical picture that emerges is that couplings that 
are irrelevant tend to deform the Fermi surface in order to become more relevant (irrelevant couplings 
being those that do not exist at vanishing excitation energy because of kinematical constraints 
attached to the Fermi surface). These insights are incorporated in a renormalization group approach, 
which allows for a simple approximate computation of Fermi surface deformation in quasi one- 
dimensional electronic conductors. We also analyze flow equations for the effective couplings and 
quasiparticle weights. For systems away from half-filling, the flows show three regimes corresponding 
to a Luttinger liquid at high energies, a Fermi liquid, and a low-energy incommensurate spin-density 
wave. At half-filling Umklapp processes allow for a Mott insulator regime where the dressed Fermi 
surface is flat, implying a confined phase with vanishing effective transverse single-particle coherence. 
The boundary between the confined and Fermi liquid phases is found to occur for a bare transverse 
hopping amplitude of the order of the Mott charge gap of a single chain. 
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I. INTRODUCTION 

One of the striking results obtained in the last decade 
on strongly correlated electronic systems is the coexis- 
tence of a notion of Fermi surface and of strong devi- 
ations from the predictions of Fermi liquid theory for 
many low-energy properties. This has been extensively 
studied experimentally for high-temperature supercon- 
ducting cuprates, where angular resolved photoemis- 
sion spectroscopy (ARPES) has revealed the presence 
of Fermi surface arcs, even in the underdoped regime 
which is characterized by the pseudo-gap seen with most 
low-energy probes.- Although these systems exhibit in- 
termediate or even strong electron interactions, they 
have triggered many theoretical works using perturba- 
tive toolsiSi^ 

At the beginning of any perturbative analysis, the 
shape of the Fermi surface is crucial in determining which 
couplings survive in an effective low-energy description."" 
For most crystalline materials the absence of continu- 
ous rotational invariance allows for a deformation of the 
Fermi surface away from the bare free electron Fermi sur- 
face, as interactions are switched on. In many metallic 
systems this effect is not expected to play much role be- 
yond usual renormalizations of effective parameters of 
band theory. But in some situations, like the vicinity of 
a Van-Hove singularity, the presence of a nesting vector, 
or for strongly anisotropic conductors, it seems essential 
to understand how to compute the dressed Fermi surface, 
since it is the relevant object for the construction of an 
effective low-energy theory. 

In the case of quasi one-dimensional (quasi ID) sys- 



tems, this Fermi surface deformation is intimately con- 
nected to the widely studied notion of transverse co- 
herence. Experimental and theoretical investigations 
converge towards a description in terms of almost un- 
coupled Luttinger liquids along the chains, at high 
enough energies. ^'^ At low energies, optical conductivity 
measurements^ have shown the existence of two types of 
behaviors: either the system remains confined in a Mott- 
insulator phase (in the TMTTF compounds) or the trans- 
verse hopping of electrons takes over and establishes a 
long-ranged transverse phase coherence, leading to a two- 
dimensional (2D) Fermi liquid phase (for the TMTSF). In 
the latter case the dressed Fermi surface remains warped 
while in the former it becomes completely flat under the 
effect of sufficiently strong interactionsAiS 

Because of their difficulty, precise computations of 
Fermi surface deformations for model systems have been 
undertaken only recently. A direct numerical evaluation 
of the electron propagator to second order in interaction 
has been performed for the 2D Hubbard modeliiiti^ Sim- 
ilar studies have also been carried for more phenomeno- 
logical models where electrons are scattered by dynam- 
ical spin fluctuationsii2*i^ Although these computations 
yield valuable physical understanding of the processes in- 
volved in the Fermi surface deformation, they suffer from 
at least two serious problems. First, they identify the 
dressed Fermi surface with the locus of points in fc-space 
for which the dressed quasiparticle energy is equal to the 
(interacting) chemical potential, which is of course cor- 
rect. But this does not imply that the imaginary part 
of the self-energy vanishes on this surface and for fre- 
quencies equal to the chemical potential. Therefore this 
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procedure does not lead to a picture of asymptotically 
stable quasiparticles at low energies. This remark is valid 
in the zero temperature approach, which is the only one 
we are using in this paper, because of its conceptual sim- 
plicity. Second, this problem is not cured while going 
to higher orders in perturbation theory. Furthermore, 
some new problems arise (namely infrared divergences) 
at these higher orders for both zero and finite tempera- 
ture formalisms. 

The underlying assumption of the standard perturba- 
tion scheme as used above is that one can generate the in- 
teracting ground-state by adiabatically switching on the 
interactions, starting from the non-interacting ground- 
state. This has to be questioned for large systems for 
which the ground-state lies at the edge of an energy con- 
tinuum. Because of this, the perturbation algorithm act- 
ing on various excited states of the original systems, as- 
sociated to different shapes of the Fermi surface, has the 
possibility to generate energy levels' crossings. This im- 
plies that the seed state to be used in perturbation the- 
ory is not known a priori, when interactions do deform 
the Fermi surface. This difficulty has been pointed out in 
the sixties by Kohn and Luttingerfi^ and also Nozieresii^ 
These ideas have been revived recently in a mathemat- 
ically rigorous frameworkJi The conclusion of all these 
works is that a sound formalism is obtained when one 
works with a bare propagator which singularities are 
pinned to the dressed Fermi surface. This is achieved 
in practice by the introduction of counterterms, which 
have to be computed order by order in perturbation the- 
ory. The main difficulty in practical implementations of 
this philosophy (which may be called renormalized per- 
turbation theory) is that it provides only an implicit de- 
termination of the dressed Fermi surface, since this algo- 
rithm expresses the bare Fermi surface as a function of 
the dressed one. Although formally this connection has 
been proved to be invertiblcfi^ this remains a formidable 
task which has never been, to our knowledge, practically 
undertaken. Note that the necessity to use these coun- 
terterms is not a pathology of the zero temperature ap- 
proach. It also appears in the Matsubara formalism at 
finite temperature which is the one used in the rigorous 
works just described. 

As a first step towards the realization of this program, 
several groups have performed self-consistent computa- 
tions. Their basic principle is to start with a trial Fermi 
surface, which is adjusted so that it matches with the 
calculated Fermi surface. A first example follows directly 
the standard Hartree-Fock method.^- It has been applied 
to the 2D Hubbard model in the presence of second- 
neighbor hopping and nearest neighbor interaction, and 
the possibility of a change in Fermi surface topology 
(from hole-like to electron-like) has been observed. A 
rather sophisticated scheme has also been developed by 
Nojiri^ in which the self-energy is self-consistent ly com- 
puted from the corresponding second order Feynman di- 
agram. This work addressed the simplest 2D Hubbard 
model with on-site interaction for which the Fermi sur- 



face deformation was found to be very small and to pre- 
serve the Fermi surface topology. Note that the quanti- 
tative difference between this self-consistent scheme and 
a standard perturbation theorjiiiiiS appears to be small. 

In spite of their merits, these approaches lack the abil- 
ity to keep track of the growth of some effective couplings, 
as the typical energy scale is lowered. These effects play a 
crucial role for the 2D Hubbard model near half-filling, or 
for quasi ID conductors. A natural way of handling these 
trends is to use a renormalization group (RG) approach. 
Several groups have incorporated the RG methodology 
in the computation of the dressed Fermi surface i^i^iiSiSi 
Similar studies have also been carried for two coupled 
chains where the Fermi surface reduces to four Fermi 
pointsi ^^i^^'^"^ Our understanding of these works is that 
they always begin with a known bare Fermi surface and 
compute the evolution of the effective Fermi surface, as 
the high-energy cut-off is gradually decreased. Although 
this is very reasonable on physical grounds, we may won- 
der whether this fits with the general rigorous analysis 
described in the last but one paragraph. We believe there 
are two ways to combine the corresponding requirements 
with a RG approach. The first one uses the renormalized 
perturbation theory described above, with a running en- 
ergy cut-off. After the usual mode integration in a small 
energy shell, the kinetic term in the effective action is cor- 
rected to preserve the shape of the dressed Fermi surface. 
In the process of integrating the RG flow, one has to keep 
track of and sum all these counterterms to obtain the bare 
Fermi surface as a function of the dressed one. Alterna- 
tively, one would fix the bare high-energy theory, and 
perform the mode integration is such a way that modes 
being integrated out always remain at a finite distance 
from the fiowing Fermi surface. But then one has to en- 
sure that all modes are integrated over exactly once with 
a uniform weight. This is indeed possible but requires 
some slight modifications of the Wilson-Polchinski usual 
RG equations. We believe the practical implementation 
of cither approach remains to be attempted. 

The bulk of this paper is composed of three sections. 
Sec. ^begins with a general discussion of some difficul- 
ties with the standard perturbation theory. We then de- 
velop a physical understanding of the driving force that 
deforms the Fermi surface on the basis of a simple vari- 
ational calculation for a system of two spinless chains. 
The main insight gained here is that the couplings which 
tend to deform the Fermi surface are those for which ex- 
ternal momenta of in and out going particles can not be 
simultaneously taken on the Fermi surface, because of 
momentum conservation. In the RG language, these in- 
teractions are usually called irrelevant. We finally estab- 
lish the equivalence between this procedure and a stan- 
dard renormalized perturbation theory where the dressed 
Fermi surface is fixed by counterterms. The reader inter- 
ested in more technical aspects is referred to Appendices 
IXI IbI and El (the first two begin with some simple first 
order calculations on the system of two spinless chains, 
whose results can be compared to the ones obtained in 
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Sec. In Sees. IIIll and II VI we show how the RG can 
be implemented in the study of quasi ID systems. We 
want to emphasise that we have not made use of a single 
RG scheme, but of two coupled RG schemes. We de- 
scribe our motivations for performing such a study in 
Sees, nil B II and IIV Al but let us very briefly explain 
what they are, before coming to a more detailed descrip- 
tion of Sees, mil and IIVI The field-theoretical RG in the 
spirit of Gell-Mann and Low^ is a simple but powerful 
way of computing low-energy properties of systems de- 
scribed by a renormalizable field-theory. This is why we 
adopted it for this purpose (this method is discussed in 
detail in Appendix ID)i. However, it cannot be used to 
compute the dressed Fermi surface, for the simple rea- 
son that the Fermi surface is defined as the locus of the 
zeros, in fc-space of the inverse propagator evaluated at 
zero frequency. There is thus no low-energy scale v that 
can be varied to get RG equations as is done for example 
for the low-energy vertices, when relating the values of 
these vertices at two different scales ly and . However, 
one can use the approach known under the name cut-off 
scaling, and developed by Solyom.'^^ This RG does not 
suffer from the limitation just described, because it is the 
high-energy cut-off and not the low-energy scale that is 
varied, and we have used it for the computation of the 
dressed Fermi surface. The high-energy part of the flows, 
in which the Fermi surface deformation takes place, is 
thus described by the cut-off scaling. The dressed Fermi 
surface that one obtains in this way then serves as an in- 
put parameter for the field-theoretical RG which governs 
the low-energy part of the flows. Let us say that RG flow 
equations appear neither in Sec. IIIII nor in Sec. IIVI but 
they all have been gathered in AppcndixlEl In Scc. lIIII we 
set up the cut-off scaling approach for the study of Fermi 
surface deformations in a quasi ID system of weakly cou- 
pled electronic chains. In order to make the ideas more 
concrete, this method is then applied to the simplest pos- 
sible example, and we end the section with a comparison 
to other methods that can be found in the literature. We 
then turn to numerical investigations, that are presented 
in Sec. IIVI for short range, Hubbard-like, repulsive elec- 
tron interactions. Sec. IIVBI deals with considerations 
about systems away from half-filling which exhibit an 
incommensurate nesting vector for their Fermi surface. 
The flow pattern involves a high-energy Luttinger liquid 
regime, followed by a Fermi liquid at intermediate en- 
ergy, and finally a long-range ordered spin-density wave 
(SDW) phase is the stable low-energy attractor. Special 
attention has been given to the scale and transverse size 
dependence of the quasiparticle weight. We then focus 
on the half-filled (and nearly half-filled) case in Sec. IIV Gl 
where Umklapp processes may drive the system into a 
confined low-energy phase and pin the SDW on the crys- 
tal lattice. In particular we study the cross-over between 
the confined and the Fermi liquid regimes. It is shown 
to occur for bare values of the inter-chain hopping of the 
order of the ID Mott charge gap. 
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FIG. 1: Schematic energy level pattern as a function of in- 
teraction strength for a conducting Fermi system. Different 
levels correspond to different choices for the Fermi surface of 
the non-interacting system. The left figure represents what 
happens in the standard perturbation theory, where the level 
repulsion at avoided crossings can not be resolved, so that one 
obtains a non-adiabatic evolution of the system wave-function 
as interactions are increased, therefore generating an excited 
state. On the right we represent the effect of applying the 
standard perturbation theory in a finite size system. In this 
case an adiabatic generation of the interacting ground-state 
is possible. 



II. COMPUTING THE SHAPE OF THE FERMI 
SURFACE: VARIOUS DIFFICULTIES AND 
THEIR RESOLUTION 



A. General considerations 

As emphasized in the Introduction, the computation 
of the dressed Fermi surface in an interacting metallic 
state encounters some obstacles because of the presence 
of a continuum of low-lying energy states in the immedi- 
ate vicinity of the non-interacting ground-state. This has 
been already discussed in a very inspiring paper by Kohn 
and Luttingerii^ There, they have shown that the stan- 
dard Brueckner-Goldstone perturbation theory for the 
ground-state energy is not consistent with a careful pro- 
cedure of taking the zero-temperature limit of the total 
energy computed in the grand-canonical ensemble. They 
interpret this failure in terms of the pattern of energy 
levels of an interacting Fermi system as a function of the 
interaction strength. When the shape of the Fermi sur- 
face changes, a deep reshuffling of the spectrum takes 
place, leading to a huge number of level crossings. A 
simple illustration for this is given on Fig. □ At this 
stage, it is important to distinguish between two situa- 
tions, which have both interesting physical realizations. 
For some simple models, such as a ladder of interacting 
spinless fermions, or a single chain of spin 1/2 electrons, 
the total number of particles of a given species (transverse 
momentum in the ladder case, or the z component of the 
spin for spin 1/2 electrons) may be conserved. As a re- 
sult of this symmetry, the level crossings just mentioned 
are an essential feature of the exact many-body spec- 
trum. In more general situations, these level crossings 
appear at any finite order in a perturbative computation 



4 



of the spectrum as a function of interaction strength, al- 
though they are expected to disappear in an exact treat- 
ment for a finite-size system. Let us first concentrate 
on the former case for a while, since it shows dramat- 
ically why and where difficulties arise. In such situa- 
tions, the conventional assumption often made in many- 
body computations does not hold. It states that one 
can get the interacting many-body ground-state by adia- 
batically switching on the interactions, starting from the 
non-interacting ground-state. A trivial example where 
the adiabatic switching procedure most often generates 
an excited state is provided in the case of the free Hamil- 
tonian: 

i^A = 5I^A(fc)ct(fc)c(fc), (1) 

k 

where we arbitrarily split e\{k) in two parts: 

eA(fc) -£o(fc) + A£i(fc). (2) 

This induces a decomposition of H\ as a sum H\ = 
Hq + XHi, where Hq is the "unperturbed" Hamiltonian, 
and XHi the perturbation. Since Hq and Hi commute, 
the eigenstates of Hx do not depend on the strength A 
of the perturbation. But energy levels as functions of 
A arc free to cross, so the initial ground-state (i.e. for 
A = 0) becomes in general an excited state for finite A. 
This is refiected on the computation of the single par- 
ticle Green's function in the zero temperature formal- 
ism. Starting with the "bare" propagator G*^"'' (fc, = 
uj — £o(fc) -I- i?7Sgn(eo(fc) — /io): the conventional algo- 
rithm yields a "dressed" propagator G^-^^^ {k,uj)~^ — uj — 
e\{k) A- i?7Sgn(eo(fc) — Mo) instead of the correct result: 
G(^)(fc, w)"^ =L0 - £x{k) + irisgn{e\(k) - nx), where 
and fix denote the bare and the dressed chemical poten- 
tials respectively. Note that the problem would appar- 
ently disappear in a finite temperature approach using 
the Matsubara formalism. However, Kohn and Luttinger 
have shown that special care is needed in taking the zero 
temperature limit, since they have found a class of dia- 
grams (they have called them anomalous diagrams) for 
which the zero temperature limit and the infinite vol- 
ume limit do not commute. Taking the former limit first 
yields a vanishing contribution for those diagrams, and 
therefore the wrong result of the standard zero temper- 
ature formalism is obtained. The correct result for an 
infinite system is obtained by taking the other order of 
limits, where anomalous diagrams do provide finite con- 
tributions. 

For this reason, and also given the conceptual interest 
of this problem, we shall use only the zero temperature 
formalism throughout this paper. In this framework, a 
natural way to circumvent this problem with level cross- 
ings is to start the standard perturbation algorithm with 
any arbitrary eigenstate of the non-interacting Hamilto- 
nian Hq. Intuitively, we believe in most cases it is suffi- 
cient to choose an initial state where the locus of occupied 
single particle states is singly connected (i.e. it has no iso- 
lated particle- hole excitations from the viewpoint of Hq), 
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FIG. 2: Examples of possible initial states for the pertur- 
bation algorithm. These states are Slater determinants with 
occupied single-particle states depicted by the dashed areas in 
fc-space. The dashed line denotes the non-interacting Fermi 
surface. In state (a), the Fermi surface is deformed, but no 
additional particle-hole excitations are present, unlike in state 
(b). 

but with a deformed Fermi surface, as shown on Fig. |21 
The selection of the correct initial state is performed by 
minimizing the total energy of the dressed state it gen- 
erates, after switching on the interactions. An example 
of this procedure is given below (Sec. IIIB|I for a simple 
two-chain model. 

For practical purposes, it is important to note that 
this approach may also be implemented through a per- 
turbative computation of the single particle Green's func- 
tion. Instead of using the free propagator G^^\k, u))^^ = 
w— £o(fc)+i?/sgn (eo(fc) — fJ-o), we should first make a guess 
for the dressed Fermi surface. This allows us to define a 
bmction $(fc) such that <I>(fc) = 1 if fc does not belong 
to the trial Fermi sea, and $(fc) — —1 ii k belongs to it. 
The locus of points in fc-space where $(A:) jumps from 
— 1 to +1 is our trial Fermi surface, and points on this 
set will be generically denoted as kp in the present dis- 
cussion. The corresponding bare propagator to be used 
in Feynman graph expansions is: 

G^°\k,Lu)-' ^uj-eQ(k)+iff^k). (3) 

As usual, the dressed propagator is obtained as 
G^{k,uj)^^ = Lo — £o{k) — I]$(/c, w), where the subscript 
$ in I]$(fc, Lj) is to stress the influence of the choice of a 
trial Fermi surface encoded in the function If this trial 
Fermi surface is the correct one for the interacting Fermi 
system, we expect the self-energy satisfies the following 
well-known conditions: 

i) There exists a well defined chemical potential so 
that for any kp belonging to the trial Fermi surface, we 
have: 

H-eQ{kp)-^j:^{kp,fi) = 0. (4) 

ii) The inverse life-time of "quasiparticles" vanishes on 
the trial Fermi surface so that: 

3S]<i.(fci.,M) = 0. (5) 

Of course, these conditions are not satisfied for most trial 
Fermi surfaces, as the reader will immediately notice on 
simple examples. We have checked on several examples 
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that both procedures (i.e. minimizing the total energy, 
or satisfying conditions i) and ii) on the dressed single 
particle propagator) yield the same dressed Fermi sur- 
face. In Appendix ^ we provide a formal proof of this 
equivalence, first in the finite volume case, and then in 
the case of an infinite volume. When the choice of $ is 
not the correct one, it is impossible to satisfy both con- 
ditions i) and ii) simultaneously. In the case of standard 
perturbation theory $ is taken to be ^^^^ corresponding 
to the bare Fermi surface, obtained from ffgiiiiiS The 
dressed Fermi surface is assumed to be determined from 
an equation which resembles condition i), namely: 



(6) 



But doing this yields two severe flaws: as shown in Ap- 
pendix ^ this does not generate the same dressed Fermi 
surface as the two procedures presented above and ar- 
gued to be the correct ones do. Furthermore, in pertur- 
bation theory, 3E$(o) (fc_F, /i) changes sign on the non- 
interacting Fermi surface, and for ld equal to the non- 
interacting chemical potential //q ^ shown in Appendix 
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This discussion holds clearly in the case where energy 
level crossings associated to various initial shapes of the 
Fermi surface are protected by some symmetries of the 
full Hamiltonian, as stated at the beginning of this sec- 
tion. Here, we would like to emphasize that a similar 
qualitative picture also holds in a more generic situation. 
On general grounds, we expect that energy levels of a 
finite system do not cross as the interaction strength is 
increased. This is the famous phenomenon of energy level 
repulsion which plays a key role in the field of "quantum 
chaos" (see for instance the book by Gutzwiller^-). So, 
standard perturbation theory starting from the unper- 
turbed ground-state is expected to generate the correct 
interacting ground-state for a finite system. However, to 
get the full single energy level resolution in the spectrum 
with all the avoided level crossings clearly requires going 
to very high orders in perturbation theory. Instead, in 
most many-body computations, we first get formal ex- 
pressions for various quantities such as Green's functions 
for a chosen finite order in powers of the interaction, 
and we most often take the thermodynamical limit he- 
fore summing the perturbation series. We believe this 
procedure is most likely to generate in the end an excited 
state of the interacting system, although the seed of the 
perturbation series is the non-interacting ground-state. 
This belief is confirmed by the simple computations in 
Appendix^ which do not require any special symmetry 
of the full Hamiltonian. 



B. Two chains of spinless fermions: Energy 
minimization 

1. Model and notations 

Let us first focus on the simplest possible model ex- 
hibiting the features described previously: a system of 
two chains of interacting spinless fermions. We will as- 
sume this system to be anisotropic, described by a tight- 
binding Hamiltonian, with a hopping t\\ along the chain 
much larger than the transverse hopping tj^. Hence, 
we have two bands, named by the transverse momen- 
tum they correspond to, i.e. (bonding) and tt (anti- 
bonding). We suppose the filling is such that both bands 
are partially filled. We will furthermore focus on the low- 
energy properties, so that we can linearize the spectrum 
around the four Fermi points, giving rise to four types of 
fermions: (R, 0), (R, tt), (L,0) and (L,7r). As usual, we 
extend the spectrum for arbitrary momenta. The low- 
energy free Hamiltonian is thus given by: 



EE (^) 

k /=0,Tr 

^(0)+40)(fc_40))l clAk)c^^,{k) + 



M<"'-4j(fc + 4°l) ci Ak)c^j{k) 



In the above expression, all the superscripts de- 
note free quantities. ^^^^ is the chemical potential, 
and A:p°j the Fermi velocity and momentum on chain /. 

4 i{k) is the creation operator of a right fermion on chain 
/, with parallel momentum k. The sum over k is to be 
understood as an integral for a system in the thermo- 
dynamic limit. In all that follows, we will simplify the 
problem and suppose that the Fermi velocities for both 
branches are equal, and they will simply be denoted as 

We shall also make simplifying assumptions about the 
interactions. Thus, the only low-energy interaction pro- 
cesses we will be interested in, are of the forward scatter- 
ing type ((72), classified as A, B, C, D, and F. They are 
represented on Fig. O We shall neglect the Umklapps, 
assuming the filling is not commensurate, interac- 
tions, involving four right or four left fermions, are also 
discarded, because we shall restrict ourselves to first and 
second order effects, to which these interactions give no 
contribution. In order to save space, we only give the D 
type interaction Hamiltonian: 



H.. 



(D) 



D 



E 



(8) 



k.k' 



where h.c. means the hermitic conjugate. 
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FIG. 3: Selected low-energy interactions for the two-chain 
model. 



2. First order 

We will here compute the energy to order one in the 
usual quantum mechanical perturbation theory, of eigen- 
states obtained from two types of free eigenstates. The 
first ones, denoted as |0; /cf.o, fcF,7r)o, are free states for 
which the bonding (respectively anti-bonding) band is 
filled up to fcp^o (respectively /cf.tt)- The ground-state 

of the free system is thus obviously |0; fcp^Q, fcp*'^)o. Of 
course, as the number of particles is fixed, the condition 
kF,o + kF,Tr = + k^\. must be satisfied. As we wish 
to understand what happens if one adds a particle to 
the system, we will also consider states that are simply 
obtained from the first ones by adding a particle of mo- 
mentum q on branch or tt (with q ^ k-p^o or g ^ ^F,7r)- 
We will refer to these states as \l,q, 0(7r); fcp.o, fcF,7r)o- We 
shall neither consider states with one hole, nor states with 
an arbitrary number of particles or holes. 

First of all we can compute the energies of these 
states, in the non-interacting case. Of course, be- 
cause our linearized dispersion relations have been ex- 
tended to include infinitely many single-particle states, 
there is strictly speaking an infinite particle density in 
these Dirac seas, which yields divergent expressions for 
the total energy. We will regularize these divergences 
by putting an ultra-violet cut-off Aq on the momenta, 
around the four free Fermi momenta, as shown on Fig.^ 
for one band. For the sake of simplicity, we work in the 
thermodynamic limit, and after a bit of algebra we find: 

i?(°)(0;fcF,o,A:F,.) = -(2m(°)Ao - ^^A^) 

TT 

(0) T 2 

+ ^ (fcp.o - 4"o) > and (9) 
q, 0(7r); fcF,o, kF^.) - E<-°^ (0; fep^o, fcp,.) 



„(o) 



(0) 



(10) 



It is obvious that the minimum of the energy is obtained 
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FIG. 4: Here we show how the ultra-violet cut-ofF is chosen 
around the free Fermi surface, for one band. 



for the free Fermi surface. The value of /i^"^ does not play 
a role here since we have fixed the total particle number. 

To order one in the couplings, it is well known that the 
energy of a free state is simply shifted by the mean value 
of the interaction for this state. As a consequence, the 
D and F couplings do not give any contribution. They 
will only start playing a role to second order. It is a very 
simple matter to check that: 



ASW(0;fcF,o,fcF,.) 



(2^)2 



A 



(Ao + (fcF,0 - 4°o) 



+B (^Ao - (fcF,o - 4°o)^ 
-2C (Ao + (fcF,o - 4°o)) (Ao - (fcF,0 



(11) 



r.(0)> 



1 

'2^ 



AS^i) (1, q, 0; fcF,o, A:f,.) = A^^^' (0; fcp.o, fcp,.) 
A (Ao + (fcF.o - 4'!o)) 

+C (Ao - (fcF.o - 4°o)) 
A£;(i) (1, g, tt; fcp.o, fcp,.) = A^^^'^ (0; A:f.o, ^f,.) 



(12) 



1 

2^ 



B (^Ao - (fcF.o 



(13) 



-C 



(Ao + (fcF.o - fcl.°|,)) 



We have used the conservation of the number of particles 
so that the above expressions are expressed only in terms 
of the Fermi momenta on branch 0. Thus we minimize 
the energy £;(!) = £;(0) + A£:(i) simply by requiring for 
its derivative with respect to fcp.o to vanish. This yields: 
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!,(!) 
^F,0 



l(0) 



(B-A)- 



Ao 



A + B-2C 



(0) 



(0) 



(14) 

Let us show how the chemical potential can be com- 
puted, using the energies of the states with one added 
particle. First of all, we notice that the expressions 
A£''^^-'(q) are independent of q. It implies the energy 
for adding a particle to the system on branch (tt) will 
be minimal if q is as small as possible, i.e. q — kp Q 
{q — fcp.Tr)- This confirms that fcp.o and fcF,7r sue the 
actual Fermi momenta. Now if we require this minimal 
energy to be the same on the two branches, equal to the 
renormalized chemical potential, we obtain the two fol- 
lowing conditions: 

^w = ^(") + 4°)(fcp^o-4°)) 

-|-A£;(i) (1, q = fcF,o, 0; fcF,o, kp,^) (15) 
-A£;(i)(0;fcF.o,fcF..), 

;.W=M(")+4°)(fcF,.-4°l) 

+AE^^\l, q = fcF,., tt; fcF,o, kp^.) (16) 
-A£;(i)(0;fcF.o,fcF..). 

One can check these equations give the deformation H14(l 
of the Fermi surface. This is physically desirable. Indeed, 
imposing that the minimum energies to add one particle 
on one branch or the other are identical, should be equiv- 
alent to the requirement that taking two particles at the 
Fermi surface on one branch and putting them at the 
Fermi surface on the other branch costs nothing (in the 
thermodynamical limit). Finally we find the chemical 
potential: 



,(1) - „(o) 



An 



H^") + {A + B + 2C)-^ 



An 



(17) 



-{B-AY^ 



Ao 



A + B-2C 



Anv- 



(0) 



Anv- 



(0) 



To conclude this section about first order computa- 
tions, we show two figures of what would happen for 
a total energy of the following simplified form: iJ^^^ = 

{kpfi - k^pl)^ + (A - B)(fcF,o - 4°o)- Fig. El illustrates 
the level crossings: we represent the energy as a func- 
tion of {B — A) (assumed positive), for various values of 
(fcF,o ~ ^f'o)- Fig. 13 proposes an alternative vision of the 
same thing (see the caption). 

Eq. (|14|l shows us that the deformation of the Fermi 
surface at first order is due to the difference between the 
couplings on the branches: if A — B, no deformation 
takes place. We can understand the sign of the deforma- 
tion very simply. Suppose the fermions repel each other 
(i.e. the couplings are positive), but that the repulsion is 
bigger on chain tt for example: B—A > 0. It is then natu- 
ral, in order to lower the energy of the system, that some 



0.02 



0.01 



Em 



-0.01 



-0.02 



-0.03 




0.35 



FIG. 5: Energies as functions of (fcp,o — fcp'o) for different 
values of B - A (0, 0.1, 0.2, 0.3, 0.35). The dashed curve 
gives the energy of the ground-state, and thus goes through 
the minima of all the different curves. 




0.35 



FIG. 6: 



Energies as functions of B 



ferent values of (fcp^o — kp q) (0, 0.05 



A, for states with dif- 
0.1, 0.15 et 0.2). The 
crossmgs show us that the Fermi surface will be deformed. 
The dashed curve is the envelope of all these curves. It is 
thus the energy of the interacting ground-state, as a function 
of the interaction. 



fermions of chain tt go to chain 0, so that (fcp.o — ^f'o) 
should be positive. This is indeed what we find. We will 
now see that things are different at second order: the 
D couplings tend to flatten the Fermi surface, whatever 
their sign, and without having to invoke a difference be- 
tween two couplings. 



3. Second order and further 

We shall now discuss in detail perturbation theory to 
second order (this notion of order being simply the num- 
ber of vertices in the corresponding Feynman graphs), 
and we will also see that some problems arise to third 
order and beyond, showing that another perturbation 
scheme is needed. 

We have already seen the effects of A, B and C inter- 
actions on the Fermi surface's shape to first order. We 
let the reader check that these three interactions play no 
essential role in the deformation of the Fermi surface at 
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second order. Indeed, when we compute the energy of the 
states jO; /cp.Oi fcF.7r}oi if we only keep contributions that 
diverge when Aq oo, we get a quantity that is pro- 
portional to Ag, and independent of the dressed Fermi 
momenta. Only finite terms do depend on the dressed 
Fermi momenta. We shall thus neglect these contribu- 
tions, and focus on D and F interactions. 

Let us begin with the effect of D interaction, which 
is the only one that does not exist at zero energy if the 
Fermi surface is not strictly flat. This is due to the con- 
straint of momentum conservation, and is easily visual- 
ized from Fig. |21 Second order perturbation theory tells 
us that the eigenenergy of an eigenstate |n), obtained 
from the free eigenstate |n)o, is obtained by shifting the 
free eigenenergy of a quantity: 



51 p(0) 



VnkVkn 



(0) 



where Vnk = o(n\V\k)o, 



(18) 



and where V is the interaction potential. This formula 
involves energy denominators. If these become smaller, 
the energy will decrease. When D interactions are con- 
sidered, we understand from these considerations that 
they will tend to flatten the Fermi surface, because this 
will allow for smaller energy denominators. That this is 
true can be checked by explicitly computing the energy 
shift, which is found to be: 



,(0) 



(27r) 



Ao 



«F,0 



«F,7r| 



(19) 

We stress that this result has been found computing ((TH|l . 
keeping only terms that are divergent when Aq ^ oo and 
that depend on the Fermi momenta. If only D terms 
are considered, it is now easy to show that the free and 
renormalized Fermi momenta are linked by the following 
formula: 



Afc. 



(0) 



D 



2Trv- 



(0) 



In 



Aq 

|A/CF 



, (20) 



where we have set Afcp = ^f.o ~ ^f.tt (and the same for 
free quantities) . This clearly shows the tendency towards 
the flattening of the Fermi surface, induced by D terms. 

What about F interactions? Perturbation theory at 
second order is divergent in the low-energy limit. In- 
deed |0; /cp.o, fcF,7r)o states that are not the free ground- 
states, are coupled to a continuum of excited states com- 
posed of two particles and two holes, which have kinetic 
energies arbitrarily close to the one of the seed state 
|0; /cf,o, fcF,7r)o- This yields energy denominators that are 
very small in absolute value, and even zero. In the self- 
energy formalism (constructed from an excited state, see 
Appendix IXI for details), this problem is regularized by 
the imaginary parts Hrj in the self-energy approach. For 
the minimization of energy, one can similarly define the 
divergent integrals with a principal part, and one finds 



A', n >. 



k^uj k + q^w + Wq k+q^uj + ujij k^i 



k' - q, to'' - lo', 

FIG. 7; Example of non-skeleton diagram giving rise to in- 
frared divergences. 



the same results as in this self-energy version. But this in- 
frared divergence is only the first one, and is not the most 
problematic. Things become worse and worse for higher 
orders. This has already been discussed by Feldman, 
Salmhofer and Trubowitaii so that we shall be brief. In 
the language of Feynman diagrams, the divergences come 
from repeated self-energy insertions, or to say it differ- 
ently, with non-skeleton diagrams. An example of the 
lowest order diagrams of this type (apart from the Kohn- 
Luttinger diagram we have already discussed at second 
order, and which is zero) is given in Fig. 13 The prob- 
lem with such a diagram is the following. Because of the 
inserted first order self-energy in the internal right prop- 
agator, we now have two right internal propagators. This 
gives a bad behavior of the integral over q around q = 
for u! = and fc = 0, once all other variables have been 
integrated out. It is clear that things get even worse if 
two or more such first or higher order self-energies are in- 
serted. We thus have to find a way of getting rid of these 
infrared problems that plague our perturbation theory. 
This is achieved by the use of counterterms, that we will 
expose now. 



C. The use of counterterms in the two-chain model 

1. Notations and first order calculation 

In order to simplify the notations, we will denote the 
Fermi velocity by vp instead of We will again sup- 
pose this velocity to be independent of the chain index 
and its renormalization will be neglected throughout this 
paper to simplify the discussion. The use of counterterms 
in interacting fermionic systems, for which the Fermi sur- 
face gets deformed, is quite old, and can for example be 
found in the beautiful discussion by Nozieres,-^^ where 
the reader will find more details. The main idea, which 
has been illustrated very recently(2LS£, is to take the in- 
teracting Fermi sea as the starting point of perturbation 
theory. As it is a priori unknown, we must ensure in 
the end of the calculation, that the "guessed" Fermi sur- 
face is indeed the dressed one. In order to have a good 
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k, LO 



A:, u) 



k, LO 



■O 



k, u) 



FIG. 8: Graphical representation of the chemical potential 
counterterm, at first order. 



FIG. 9: Graphical representation of the Fermi momenta coun- 
terterms, at first order. 



starting point, the most natural idea is to split the free 
Hamiltonian into two bits: one that is a modified free 
Hamiltonian with the correct interacting Fermi surface, 
and another that will be the difference between the true 
free Hamiltonian, and the modified one. We will thus 
write: 



EE 

7=0, 7r k 



[/I + v-p{k - fcpj)] c\^j{k)ci^ i{k) 



(fc) 



(21) 



with 



<cl = 5/^E4„/(fc)%,/(fc) 



7,fc 



and 



(22) 



I,k 



Hq% = "VpSki^cl^jikyj^jik) 



I,k 

-VFSkjJ2cli{k)c^,Ak). (23) 

I,k 



The counterterms Sfi and Ski are found by the require- 
ment that Hq remains the true free Hamiltonian: 



(0) 



l(0) 
^F.I 



-11 + S^Ji, with 
- k-p.i + (5fc/, with 



(5fcr — 5k 



(1) 



5k 



(2) 



(24) 
(25) 
(26) 
(27) 



Note that we have used a symbolic notation G" for the 
couplings A to F, and the sum over a is implicit. We 
stress that there is not only one counterterm for the 
chemical potential (or for the Fermi momentum of each 
chain), but an infinity, which are all the 5/i'-"-''s, for 
n = 1,2,.... The number n gives the power in the cou- 
plings of the considered counterterm. Counterterms have 
to be computed order by order, one after the other, in a 
perturbation theory. When using the counterterms, the 
Luttinger theorem simply says that: 5ki = 0, or for 

each order j: J^i^^'i^ = 0- ^ow the free (R,0) prop- 
agator with which Feynman diagrams are computed is: 



1 



U! - [^ + VF{k - kp /)] + i77Sgn(fc — fcp /) ' 

(28) 



and similarly for other types of fermions. Both real and 
imaginary parts of these propagators refer to the inter- 
acting Fermi surface. 

We shall now see how to implement the use of countert- 
erms in the perturbation theory of the two-chain model. 
For this, it is useful to associate a graphical represen- 
tation to the counterterms. This is illustrated at first 
order in Fig.|Slfor the chemical potential, and in Fig. 
for the Fermi momenta. The chemical potential countert- 
erm is represented by a square, whereas the Fermi mo- 
menta counterterms are denoted by hexagons. In both 
cases, the number written inside the symbol is the order 
n mentioned previously. Notice that for Fermi momenta, 
we do not need to explicitly write down the chain index 
/, because it would be redundant with the chain index 
of the propagators. The reader should also remark that 
counterterms for right or left fermions are exactly iden- 
tical. 

Now that the general notations have been given, let us 
see what the counterterm approach gives to first order, 
for the two chains. In all that follows, we will not use 
an ultra-violet cut-off around the free Fermi surface, but 
around the interacting Fermi surface. This will slightly 
alter the results, but it makes the computation simpler, 
without involving a qualitatively different physics. The 
tadpole diagram of Fig. ^| computed with the new free 
propagator g1°^ and the new cut-off, gives the following 
contribution to the self-energy: 



I]«(fc,^) 



eW(A:,^) = (B + G)^ 



(A + G) — , and 
27r 

Ao 



27r 



(29) 
(30) 



But we also have the counterterm contributions (dia- 



\ 



k, u) 



k, u) 



FIG. 10: First order contribution to the self-energy: the tad- 
pole graph. 
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k' — q,ui' — (Jij 



FIG. 11: Sunrise diagram contributing to the second order 
self-energy. 

grams of Figs. |H1 and : 

^^!:;io(k,^) = ^'^!Zik,u;)-s^/'\ (31) 

5^ct;R,o(fc'^) = -«F^fcr, (32) 
sJ^.(^-^) = "^F'5fc^^- (33) 

The dressed propagators G are such that they satisfy the 

Dyson equation: G^^ = Gf ^ — S — Ect- The chemical 
potential and Fermi momenta are found by requiring that 
they vanish for lo — fi and for k on the interacting Fermi 
surface, and that the Luttinger theorem is satisfied: 

GR,o(fc = fcF,o,c^ = A') = 0, (34) 
G^^,{k^kF^^,Lu^^i)^0, (35) 

J2 ^^^P = (36) 

with j = 1 here, because we're working at first order for 
the moment. It is very easy to check that one finds: 

5^(1) = _(A + B + 2C)^, (37) 
47r 

This is fully compatible with equations H14I) and (|17|l . 
except for second order terms that we do not find here, 
because we have changed the way we choose the cut-off. 



k',u)' 




FIG. 12: Kohn-Luttinger diagram contributing to the second 
order self-energy. 



2. Second order calculation with counterterms and 
next-order considerations 



As for the first order calculation, we have "usual" 
contributions to the self-energy, namely the sunrise and 
Kohn-Luttinger diagrams of Figs. ^] and ^1 We also 
have "pure" counterterms contributions, as in Figs. |S| 
and |51 with the index 1 replaced by an index 2. But 
now we also have two "mixed" contributions, involving 
counterterms of the previous order, shown in Figs. ll3| and 
n~i1 In fact, these two graphs vanish, for the same reason 
the Kohn-Luttinger graph vanishes. This is consistent, 
because there is no divergence from the Kohn-Luttinger 
graph to cancel. 



Before studying the sunrise graph, let us see how use- 
ful the counterterms are for third order graphs, on the 
example of Fig. [7| It is now obvious that it will be com- 
pletely canceled by the same graph, with the inserted 
tadpole replaced by the two first order counterterms. No- 
tice that the fourth order graph, consisting of still the 
same graph, with an inserted sunrise instead of a tadpole, 
would not be canceled by the graph with inserted second- 
order counterterms. The reason is that the sunrise is 
frequency and momentum dependent, but the countert- 
erms are not. However, the counterterms allow for the 
infrared divergence cancellation obtained at zero external 
momentum and frequency. 



The sunrise is easily computed, and one gets the fol- 
lowing contributions (the interaction index G refers to 
the interaction associated to the two black dots in the 




FIG. 13: "Mixed" contribution to the second order self- 
energy. This graph is a tadpole, with an insertion of the 
first order chemical potential counterterm. 
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k, ui 



k, w 



FIG. 14: "Mixed" contribution to the second order self- 
energy. This graph is a tadpole, with an insertion of the 
first order Fermi momentum counterterm. 



sunrise): 



Sj^ Q.(5(fc = fcp^O + K, = /i + I/) 



1 / G 



4 V 27rwF 



(v — v-pn) In 



(2^;fAo 



for G = A, C, F, 



(2) 



1 / G 



{v — vyh) In 



for G== B,C, F; 



(39) 



(40) 



- (tj:^) Afcp In (I AA:f|/Ao), which is nothing but 

The dressed (R,0) propagator (and others as well) can 
finally be deduced from all this: 



G 



(2) ~ 

Q (k ~ A;f,o + ^ = /i + i^) = — VFf^ 



(—] 



C 



X (i/ — vfk) In 

--(— V 

4 V 27rwF / 



27r?;F 

|^^^(^F^)"| 

(2wfAo; 



(—)' 



(j/ — vfk) In 



2 y (46) 

- [t;F(K + 2AfcF)]2| 



1 / D 



4 \27rwF 



(2t;FAA:F)ln 



y ~ [fF(K + 2AfcF)]^| 
(2i;FAfcF)2 



We could now define renormalized propagators, introduc- 
ing a wave function renormalization, and show how to 
implement a RG calculation of the dressed Fermi sur- 
face. In order not to be too redundant, we will do this 
for the more general case of N chains of spin 1/2 elec- 
trons, which is anyway physically motivated by the case 
of quasi ID systems. 



III. CUT-OFF SCALING RG CALCULATION 
FOR A SYSTEM OF TV CHAINS OF SPIN 1/2 
ELECTRONS: FORMALISM 



(2") 



1 f D 



4 V 27rwF 



X In 



[v -VY{n + 2 Afcp)] 
(2i;fAo)2 



l'2'\ 

^R7r;Z5(^ = ^F,7r + K.,UJ = il + v) 



1 / D V 
^-{- I [iy-VF{K~2AkF)] 



4 \27rwF / 
X In 



(41) 



(42) 



The second-order conditions ensuring that the trial Fermi 
surface is indeed the interacting one read: 



AA:f 



1 f D Y 
— (-2wFAfcF)ln , . 

-S^i^^^ + vpSk^^'^ = 0, (43) 

-(5/Lt(2) + vpSk'^'' = 0, (44) 
Sk^^^ + 5k^^^ = 0. (45) 



These equations lead to (5/^^^^ = 0, and to Jfc, 



(2) 



A. Setting of the model 

The free Hamiltonian is much like the one of Eq. 0, 
except that there are now N chains instead of 2, and that 
the fermions carry a spin index a: 

N 



k 1=1 a=], I 



(0) i „(o) 



„(0) 



(0). 



ik)c^jjk) (47) 



(kWiJk) 



We will in fact assume, as we did previously, that the 
Fermi velocity is independent of the chain index I, and 
that it remains unrenormalized. We will thus simply use 
the notation vp. 

As in the two-chain model, we select low-energy in- 
teraction processes. Those are of two types. The first 
one denoted by G, generalizes the interactions A to 
of the two-chain model. They are forward or backward 
scattering interactions. We shall only be interested in in- 
teractions that are invariant under spin rotations. Thus, 
we will use the charge and spin couplings G^ and G^. 
We refer the reader to our previous paper~ for more de- 
tails about this parametrization. There is however one 
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FIG. 15: Graphical representation of the spin interaction 
G%{I,J). 




J J 



FIG. 16: Graphical representation of the Umklapp interaction 



major dilFerence between the situation described in this 
article, and the one we are interested in here. Because of 
periodic boundary conditions in the transverse direction, 
all indices /, J and 6 are defined modulo the number of 
chains iV. This was not the case in our previous article, 
where the chains were obtained after considering patches 
on a nearly square Fermi surface, thus the A'^ chains had 
boundaries, and as a consequence the chains were not all 
equivalent. 

Furthermore, if the filling is not too far from one half, 
we have to consider Umklapp scatterings. These will be 
denoted by U . It is easy to convince oneself that due to 
the Pauli principle, there is no need to consider exchange 
couplings for the Umklapps. The interaction Hamilto- 
nian is thus: 



Hint ^ Hif + H^"J , with: 



(48) 



AG) ^ 27r-i;F 
i"t ~ NL 



H-i'^^T. E E E 



G'sil, J)Ir.r'lp,p' + GUI, JWr.r' ' Tp.p' 



(49) 



and 

,U) _ TTWF 

" NL ^ ^ 2^ 2^ 

I ^J.S k,k' ,q t^t' p^p' 
\^s{I,J)lr,r'lp,p' (50) 

+ h. c. |. 

The factors \/N are required to yield a good ther- 
modynamical limit. The terms have been factor- 



ized, so that the couplings are dimensionless, and this 
will suppress many 2'itvy denominators in the following. 
In the case of G couplings, the left-right symmetry re- 
quires Gf''\l,J) = G!_^^-^(J,/), and the hermiticity of 

i/int yields Gf\l, J) G"J5^(/ + S,J -5). The first of 
these relations, i.e. Us{I,J) = U^s{J,I), naturally holds 
for the Umklapps because of the Pauli principle, so that 
the interaction that destroys two left fcrmions on chains 
/ and J, and creates two right fermions on chains 1 + 6 
and J — S is present twice. The difference of a 1/2 factor 
between the Umklapps and the G interactions, is here to 
compensate this. We let the reader check that in the case 
of the Hubbard model with an interaction Hamiltonian 
^ Si "•j.T'^i.i' '-"^^ (^P ^'^ 27ruF factors) G"^ = U/2, 
G^ — ~U/2 and U =U. Because of this last equality, we 
will simply give the value of U when referring to Hub- 
bard couplings. Of course the Hubbard model, in terms 
of right and left fermions, also contains interactions, 
but these have been set to zero, for the reasons already 
given in Sec. Ill B II 

In order to make our notations for the interactions a bit 
more concrete, we show two Feynman graphs in Figs. lT^ 
and El associated respectively with G^ and U terms. 
The representation for G"^ is the same as the one for C, 
except it involves I matrices instead of a matrices. Notice 
we do not use the single dot notation as we did previously, 
because it is not suited for the Umklapps, but we have 
adopted the wiggly line instead. In these graphs we also 
show which external legs are numbered 1, 2, 3 and 4. 



B. Cut-off scaling calculation of the Fermi surface 

1. General considerations 

One of the main conclusions of Sec. ^1 is the neces- 
sity to use a renormalized perturbation theory in situ- 
ations where the Fermi surface changes as a function of 
interaction strength. In the standard many-body formal- 
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ism, this is achieved by the introduction of counterterms 
which pin the dressed Fermi surface. The outcome is a 
precise connection such as Eq. H2()(l between the bare and 
dressed Fermi surfaces which may in principle be com- 
puted to any order in perturbation theory. As noticed 
already long ago by Gell-Mann and Low,'^° it is possible 
to sum infinite classes of contributions using a renormal- 
ization group procedure. This idea has played a crucial 
role in building a consistent physical picture of quasi ID 
conductors for instancei^ 

The most general and flexible way to implement a 
renormalization approach is based on Wilson's idea of 
gradual mode elimination. Several groups have recently 
implemented Wilson's approach to the RG, expressed via 
the Polchinski equation^S*^ or its one-particle irreducible 
version>^£2, Although these equations are exact, they are 
quite complicated, since effective interactions involving 
an arbitrary number of particles are generated along the 
RG flow. Any numerical computation requires therefore 
drastic truncations in the effective action. For this rea- 
son, we have prefered to use a simplifled version of RG 
which is known as "cut-off" scaling. This procedure has 
been initiated in the pioneering work by Anderson et 
al^ for the Kondo problem, and put in a more mathe- 
matical form by Abrikosov and Migdal?"* and Fowler and 
Zawadowski^ A very extensive review on this method 
has been written by Solyom^ 

This scheme amounts to constructing a one parame- 
ter family of "bare" Hamiltonians. These are deflned 
on the single particle states whose momentum lies in 
a strip of width A away from the Fermi surface. It is 
therefore natural to parametrize these Hamiltonians as a 
function of A. Note that by contrast to Wilson's effective 
action, which includes all the possible types of interac- 
tions (relevant, marginal and irrelevant ones), the cut-off 
scaling procedure only considers relevant and marginal 
couplings. So unlike what is achieved in Wilson's RG, 
it is no longer possible to preserve invariance of the full 
set of low-energy correlation functions as A is gradually 
decreased. The cut-off scaling approach only allows to 
preserve a restricted set of low-energy observables, for 
instance the first derivatives of the two-point function 
with respect to external momentum and frequency, and 
the value of the four-point function for external legs taken 
on-shell at the Fermi level. 

Actual computations within this scheme encounter a 
new difficulty when the Fermi surface is sensitive to 
the strength of interactions. As explained in detail in 
Sec. Ill Gl the bare propagators used in Feynman graphs 
are required to be singular on the dressed Fermi surface. 
But clearly, this is not known until the whole compu- 
tation has been performed. This calls for an iterative 
procedure. For a given microscopic model (defined with 
an initial cut-off Aq), the one-particle part of the corre- 
sponding Hamiltonian defines a Fermi surface which will 
be called the Ao-Fermi surface. This is a natural first 
choice for a trial dressed Fermi surface in the iterative 
computation. One can then construct the fiow equations 



for the running bare Fermi surface (the A-Fermi surface) 
and the effective couplings using cut-off scaling. In the 
limit where A goes to zero (or at least to its minimal 
value before a phase transition occurs), the A-Fermi sur- 
face goes towards a new dressed Fermi surface, which is 
to be used as the new trial dressed Fermi surface in the 
next step of the iteration. 

On physical grounds, such a computation is expected 
to converge although we have not embarked yet in check- 
ing this statement. Instead we have tried to bypass the 
intrinsic complication of an iterative procedure by ap- 
pealing to the physical insights gained in Sec. Ill Bl de- 
voted to the energetic approach. The main ideas are the 
following: firstly, in the high-energy regime, and in the 
logarithmic approximation, the dressed Fermi momenta 
do not appear in the flow couplings' equations; secondly, 
in the low-energy part of the flow, the Fermi surface 
should not move too much, because the couplings that 
deform it are irrelevant in this regime. The first point 
will be checked on the RG equations. The second one has 
already been checked in the two-chain model, where only 
the D coupling, that does not exist at very low energies 
because of the non-flatness of the Fermi surface, deforms 
the Fermi surface. We will furthermore check it remains 
true in the case of N chains. Assuming what happens in 
the intermediate energy regime (defined by the curvature 
of the Fermi surface) is not essential, the computation 
of the dressed Fermi surface is now possible in a single 
step. Indeed, we do not need a priori knowledge of the 
dressed Fermi surface anymore, since it disappears from 
the RG couplings' flow equations in the logarithmic ap- 
proximation of the high-energy regime. The flow of the 
A-Fermi surface, will then be stopped when the cut-off 
A becomes comparable to the maximal momentum scale 
defined by the running A-Fermi surface. Note that this 
is the least controlled step of this approximated scheme, 
because the Fermi surface does not define one single mo- 
mentum scale, but rather a continuum of scales. This will 
for example prevent us from using this scheme when the 
Umklapp couplings are taken into account, in a system 
too far from half-filling. 

The couplings' fiow equations in the cut-off scaling 
are given in Appendix ^ since their derivation is stan- 
dard. We shall now focus on the RG computation of 
the Fermi surface. The basic equation is Eq. (|D9|I . The 
self-energy that appears in this equation can be found 
from Eq. (|E13I) . where, since we work in the cut-off scal- 
ing scheme, the cut-off Aq should now be replaced by 
the running cut-off A, and where the couplings are to 
be understood as running bare couplings. Given that 
Ski = kf\A) - ki, Eq. jnil) allows us to express kf\h.) 
as a function of A, the running bare couplings and the 
dressed Fermi momenta fc/. But it is the latter who are 
fixed independently of the value of A, so that it is more 
convenient to invert this relation, working only at a sec- 
ond order accuracy, to get kj as a function of A, the 
running bare couplings and the running Fermi momenta 
fc|°^(A). Asking for the invariance of kj as A is changed 
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yields the Fermi surface flow equation that we give in 
Appendix El Eq. JETSJ. 

Only the couplings for which the curvature of the 
Fermi surface is felt, i.e. for which Aka{I,J) ^ or 
Afc^(/, J) ^ contribute to the flow of kf\A). This 
is what we had already noticed in the two-chain system, 
where the D coupling was the only one to give a defor- 
mation of the Fermi surface. This confirms that only 
couplings that will be irrelevant in the low-energy regime 
contribute to the deformation of the Fermi surface. 

Finally, notice that we have not taken the first order 
contribution into account. This is not justified in the gen- 
eral case, but for the initial condition we are interested 
in, i.e. with all charge, spin and Umklapp couplings equal 
to Gg, Gg and Ub, the first order contribution vanishes. 
We emphasize this is true only in the high-energy regime, 
where the couplings will have a purely one-dimensional 
(ID) flow (because of the logarithmic approximation), 
and thus will remain equal (with one value for each of the 
three types of couplings). Once the low-energy regime 
is reached the couplings become different, so that one 
should take the first-order deformation into account. 

The self-energy at one loop is given by the contribution 
of the tadpole diagram. It is easy to see that out of the 
three couplings G'^, and U, only the G'^ couplings 
contribute. Compared to the spinless case, there will be 
a factor of 2, because of the two possible spin states of 
the propagator in the loop. We let the reader check that: 



^R,7(^ ~ ^F.-f + 1^, ^ — fJ- + J^) 



(2i;fA), (51) 



and that the corresponding first-order Fermi momenta 
counterterms read: 



(1) 



^(EGS(/,^)-^EG§(/,j)j 



(52) 



However, because we are interested in systems for 
which t± is small, i.e. for systems that are nearly ID, we 
know that all the chains are nearly equivalent, so that 
the RHS of the previous equation will be nearly inde- 
pendent of /, and thus, very small. That this is true 
can be checked on Fig. that will be described later, 
and on which the couplings Gq(/, J) are represented af- 
ter running the flow, into the low-energy regime. It is 
clear on this figure that the term J^j Gq(/, J) is nearly 
independent of /. 



2. Analytical study of the simplest example 



Let us illustrate all this on the simplest possible case, 
for which Gg = and /7b = 0. This physically corre- 
sponds to a system away from half filling, so that it is 
justified to neglect the Umklapps. Furthermore we have 



set the spin couplings to zero, which corresponds to the 
Luttinger liquid fixed point. This simplifies the flow, be- 
cause the charge couplings then remain constant along 
it. Furthermore, this will allow us to compare our results 
to those obtained in the literature, starting from decou- 
pled Luttinger liquids, coupled by a hopping term. It is 
easy to check that the general flow equation of the Fermi 
surface ()E15p can be simplified in: 
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If we denote hy k = {J2ikFj)/N = (E/ 4°!)/^ t^e 
mean value of the Fermi momenta, and if we write 
kp\{t) — k + Skp\{t), the differential equation is eas- 
ily solved and the solution is: 



<54"}(A) =<54^^(Ao) 



(0), 



(54) 



The question that now arises, is how to determine at 
what scale A* the flow should be stopped. This scale can- 
not be determined precisely in the cut-off scaling scheme, 
which is only a very simple version of the Wilsonian ap- 
proach. Only the latter approach could precisely describe 
the transition between the two regimes, which would not 
even occur at a single scale (because of the large num- 
ber of different scales K appearing in the RHS of the 
RG flow equations). We will thus adopt the simple and 
pragmatic following point of view: the flow of the Fermi 
surface will be stopped, when the biggest of these scales 
is reached, i.e. when the scale given by the difference be- 
tween the biggest and the smallest Fermi momenta (de- 
noted by Afcp''''(A)) is reached: 



A* = Afc]?^'^(A*). 



(55) 



Notice that if the Fermi surface flattens more quickly 
than the RG time decreases, this scale will never be 
reached. 

According to Eq. (|54|l . the differences between the 
Fermi momenta and the mean value are all multiplied 
by the same factor. The biggest (respectively small- 
est) momentum will remain the biggest (respectively 
smallest) along the flow. We thus have Afcp''''(A) = 
Afcj?^'^(Ao)(A/Ao)"^B'. We will stop the flow at the 
scale A* such that A* ~ AA:}?^'^(Ao)(A*/Ao)'^b'. Finally 
we get the following link between high-energy and low- 
energy momenta: 



Skpj ~ <54°1(Ao) 



Akf^''{AQ)\^ 



(56) 



A.p-(Ao)(^^l^)"^.(57) 
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As Afc™'^''(Ao) = 2tj_/t||, we obtain the result already 
found in the literature&iS 



(58) 



where a is the single-particle Green's function's ex- 
ponent: a = [Kp + l/i^p)/4 - 1/2, with Kp = 
^{l-2G%)/{l + 2G%). Perturbatively, a = G| ^ so 
that our result is indeed the same as Eq. (|58|) . to lowest 
order. Let us mention that the power-law behavior of 
Eq. H58I) has been confirmed numerically, for a two-chain 
system, using exact diagonalization techniques 

Notice that according to Eq. (|57|) . if Gg^ ^ 1, the effec- 
tive transverse hopping vanishes, and the dressed Fermi 
surface is flat. Although this result is confirmed by the 
non-perturbative (in the coupling) result Eq. H58|) after 
replacing Gg^ by a, we shall not use the perturbative 
RG in such situations that lay outside the validity range 
of this approach. 



3. Comparison to other previous results of the literature 

Before turning to numerical calculations, we shall com- 
pare our equations describing the deformation of the 
Fermi surface to some more results of the recent liter- 
ature. Let us begin with the article by Kishine and 
YonemitsufSi which treats exactly the same problem as 
ours. We shall compare our equations (|E15|) to their flow 
equation for the effective transverse hopping (Eq. (4)). 
Note that they obtained this equation from the pre- 
vious works by Bourbonnais and co-workerSiMLSS, and 
Kimurai2S (see also the review article by Firsov, Prigodin 
and Seidel"). We however choose the paper by Kishine 
and Yonemitsu because their formalism is the closest to 
ours. 

The comparison is easily achieved in two steps: first 
notice that our Fermi momenta do not flow when no in- 
teractions are turned on (which is desirable, since the 
Fermi surface should not get deformed in this case), while 
their t± flows in this case, because of the first term of 
their equation coming from the rescaling they have per- 
formed, so that we should simply forget this term if we 
want to compare our results. Second, we have to take 
the particular set of couplings they have chosen, namely 
local couplings. This is obtained when setting all charge 
couplings to the same value, and doing the same for spin 
couplings and Umklapps. Repeating exactly what we 
have done in the previous section, we find: 
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This in particular means that all Fermi momenta 5k'^\ 
will be scaled by the same factor. Next, we have to link 
our charge and spin couplings to the g-ology notation. It 
is easily checked that one simply has: G'^ = (72 — 51/2 and 



G^ = —51/2 (here all g\\ and gj_ couplings of the g-ology 
are equal because we have restricted ourselves to spin- 
rotation invariant couplings). We also have the trivial 
identification U — g^. The difference in the numerical 
factor 4 simply comes from a different normalization of 
the dimensionless couplings (we divided the couplings by 
2t:vf and they divided them by nv-p). Finally, we have 
^^max _ 2(f_L/fj[) which shows the equivalence between 
the two approaches. 

We want to stress that this equivalence relies on our 
simple approximation according to which the Fermi sur- 
face is deformed only in the high energy regime, since this 
deformation is driven by irrelevant couplings. It would 
be possible to go beyond this approximation by imple- 
menting the iterative procedure outlined in Sec. IIII B II 
To estimate the residual deformation of the Fermi surface 
induced by these irrelevant couplings in the low-energy 
regime remains an interesting open question, that could 
be addressed within the general framework discussed in 
this paper. Furthermore, such a calculation would enable 
us to take into account the deformation of the Fermi sur- 
face induced by the Hartree terms which are effective only 
when the forward scattering amplitudes significantly vary 
along the Fermi line. Such a dependence is only gener- 
ated when the running cut-off becomes comparable to or 
smaller than the natural scale associated to the trans- 
verse dispersion. 

For the sake of completeness, we shall give a simple 
and quick derivation of the RG equation, which empha- 
sizes the role of the ID chains (see note 31 of Ref.'^^), 
and explains why the exponent obtained in Eq. (|58|l is 
the ID propagator's exponent. The idea is to assume 
that the full propagator at scale A can be obtained by 
taking the corresponding purely ID propagator at scale 
A, and correcting it with the dispersion relation induced 
by the bare tr- G'^i^) = Gioi^) + 2tj_cos{k±). (In 
other words, this amounts to assume that when comput- 
ing the effective action at scale A, one puts the inter- 
chain hopping aside, so that the flow is purely ID, and 
the (unrenormalized) inter-chain hopping is reintroduced 
in the effective action at the end of the computation). 
But we can write G^^ = Zjjj(A)[[x) — eA(fc||)]- Note 
that in the previous two formulas, we have denoted by 
k± and fcy the transverse and longitudinal momenta. 
ZiD is the ID wave-function renormalization, and eA(A;||) 
is the renormalized ID dispersion relation. We thus 
get G-\A) = Z-j^iA)[uj - EAikii) + 2Zu,iA)t± cos{k±)], 
showing that the effective inter-chain hopping at scale A 
reads: t±{A) = ZiY){A)t±. The effective t± at two dif- 
ferent scales are thus proportionally related by the ID 
Z function, whose flow equation can easily be deduced 
from Eq. (jE14|) specialized to the ID case. This yields 
the correct flow equation for tj_ . Let us also mention that 
this way of taking into account the inter-chain tunneling 
has been recently adopted by Essler and Tsvelik,^"'^ ex- 
cept that they use the exact ID Green's function instead 
of the result of a perturbative RG computation. 

Let us now compare our results with those of 
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Fabrizio,^^ whose work is devoted to the two-chain model 
without longitudinal Umklapps, but with Fermi veloc- 
ity renormalization. Fabrizio used a Wilsonian RG (at 
two loops) for the calculation of the deformation of the 
Fermi surface. As one can expect, this formalism allows 
to cross the energy scales coming from the non-flatness 
of the Fermi surface (see how the flows are defined piece- 
wise in his appendix A, and for which the various RHS 
never diverge) . Our equations coincide with those of Fab- 
rizio in the high-energy regime (when his function C2 is 
expanded to lowest order in h, the dimensionless Afcp), 
and in the low-energy regime where the flow of the Fermi 
surface vanishes. The intermediate regime is of course 
different. Note for the comparison that Fabrizio's cou- 
pling gb is our coupling D. 

Finally, we would like to note that our results at one 
loop are consistent with the article of Louis, Alvarez and 
Gros,^ (see their Eqs. (9) and (13)) if we specialize these 
to the case of uniform Fermi velocities. 
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FIG. 17: Flow of the normalized charge (top) and spin (bot- 
tom) couplings, as a function of the "good" time s, for 
A'' = 8. Initially the true couplings (not normalized) are 
G" = 0.3 = — G% and the bare hopping is t±/t\\ = 0.1. We 
have indicated the value s* corresponding to the value A* (see 
Eq. rai. 



IV. COUPLED CUT-OFF SCALING AND 
FIELD-THEORETICAL RG CALCULATIONS 
FOR A SYSTEM OF N CHAINS OF SPIN 1/2 
ELECTRONS: NUMERICAL RESULTS 

A. Motivation of the use of two RG schemes 

Before we present our results for specific models, we 
wish to emphasize that one of our motivations besides 
the Fermi surface deformation was to describe precisely 
the connection between the essentially ID high-energy 
regime, and the 2D low-energy physics where the system 
is sensitive to the warping of the Fermi surface. This can 
be viewed as a complement to the RG analysis of Lin 
et ali^ who have focused exclusively on the low-energy 
side where the cut-off is much smaller than the scale as- 
sociated to the transverse dispersion. In this work, they 
could relate the high and low energy regimes without 
actually solving RG flow equations for the former since 
they assumed very weak bare couplings (so that these 
couplings were barely renormalized in the high energy 
part of the flow) . Although working with the full Wilso- 
nian effective action allows one to get through such in- 
termediate energy scales, it is not clear to us that this 
can be achieved in a reliable way with the cut-off scaling. 
Indeed, our view of this procedure is that it provides a 
simple approximation of the full Wilsonian RG, which is 
certainly well controlled when the running cut-off is much 
larger than the intrinsic low-energy scales of the system's 
dynamics. Because of this we have decided to study the 
low-energy part of the flow in the field-theoretical frame- 
work. This latter scheme heavily relies on the existence 
of an infinite cut-off limit (continuum limit) , or in other 
words the corresponding theory of ID fermions with lin- 
ear dispersion and point-like interactions is renormaliz- 
able. This statement is independent of the existence of 
intrinsic low-energy scales such as a mass term, or varia- 



tions in Fermi wave vectors with the chain index. In this 
context RG equations are obtained by relating physical 
properties measured at different running energy scales. 
To avoid confusion with cut-off scaling this running scale 
has been denoted by ly in Appendix ]n\ which presents 
some details on this approach. Since we have used a 
logarithmic approximation, the high-energy flows of the 
couplings in both cut-off scaling and field-theoretical RG 
are identical. It is an interesting question whether the 
two schemes give the same physical low-energy results or 
not. We plan to study this in more detail in a forthcom- 
ing work. 



B. A first numerical study: incommensurate 
nesting 

We will now show what information can be deduced 
from numerical computations. For this we choose to focus 
on a simple example, where the Umklapps are set to zero, 
but we still assume a perfect Fermi surface nesting. This 
situation is realized in several interesting systems as for 
instance in two dimensional molybdenum and tungsten 
bronzes. For a review, see for example the paper by Foury 
and Pouget.^'* 

We have chosen an initial condition for which all charge 
couplings are equal, and all spin couplings too, with 
G"= = 0.3 = -G". The bare hopping is t^/tn = 0.1. The 
couplings are quite large so that the deformation of the 
Fermi surface will be visible. The Fermi surface could be 
deduced analytically, but we have computed it numeri- 
cally as all other quantities. All the results are contained 
in Figs. El to 1221 The first three (respectively last three) 
of these figures have been computed with N = 8 (respec- 
tively TV = 32). The reasons for these choices are that we 
could not represent all the couplings (the first of the six 
figures) for a too high value of TV, because the number 
of couplings grows like iV'^ . But this was no problem for 
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the Fermi surface and the quasiparticle weights, except 
for a longer computation time. 

In these figures, that we shall comment one after the 
other, we have made use of some notions such as the 
norm of the couplings, the normalized couplings, and the 
adapted RG time s. All these notions, and some oth- 
ers (such as fixed directions, etc.) have been dealt with 
extensively in our previous paper^SS so we shall simply 
give the few basic definitions. The norm is the Euclidean 
norm of the coupling vector, and the normalized cou- 
plings are the usual couplings divided by the norm (we 
give an explicit formula for the charge couplings only): 



^ = J E + ^sil^ Jf + Us{I, J)\ (60) 



I,J.S 



(61) 



In all the figures in which flows will be represented (such 
as Fig. I17|l . we adopt the following convention: in the 
caption, we give the initial condition for wswa/ couplings, 
while in the figures themselves we draw the normalized 
couplings, and suppress the " in the y-axis legends. The 
reader should not get confused by this abuse of notation. 
The time s that we have used in the numerical simulation 
is defined by ds = N{t)At, and is the time adapted for 
zooming on the flow singularities. 

The first of the six figures. Fig. I17L represents the 
"field-theoretical" RG flow of the normalized charge and 
spin couplings, as functions of the RG time s. This flow 
is divided into three regions. In the first one (0 s 
s* ~ 15), corresponding to the high-energy regime, all 
charge couplings and all spin couplings remain equal. In- 
deed, in this regime, the curvature of the Fermi surface 
is not felt at all, in the logarithmic approximation we 
use. All chains are thus identical (remember we use peri- 
odic boundary conditions in the transverse direction), the 
system is purely one-dimensional, so that the symmetry 
between the chains cannot be broken. 

In this high-energy regime, the cut-off scaling flow is 
exactly the same as the "ficld-thcorctical" one, so that 
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FIG. 19: Link between the "good" RG time a and the "true" 
RG time t, corresponding to Fig. 1171 We have indicated the 
values of s* and of t* — ln(Ao/A*) corresponding to A* 



we could use the latter for the computation of the de- 
formation of the Fermi surface. After s ~ 15, the flow 
of the Fermi surface is stopped, because the scale A* (as 
previously defined by A* = Akp'^'^{A* j) is reached. The 
dressed Fermi surface is the one obtained at that scale, 
and is then used in the "field-theoretical" flow of the 
couplings. This dressed Fermi surface, and the bare one, 
are represented on Fig. 1201 As discussed in Sec. IIII B 31 
within the approximation we use, the dressed Fermi sur- 
face is still given by t± cos{k±) but with the dressed value 
of the interchain hopping. Higher harmonics for Afcp as 
a function of k±, corresponding to longer range trans- 
verse hoppings, are expected to be generated only in the 
low-energy part of the Fermi surface flow. Indeed it is 
only in this regime that effective couplings acquire a dra- 
matic dependence with respect to transverse momenta. 
But this goes beyond the scope of the simple (i.e. non it- 
erative) procedure used for the numerical computations 
presented here. 

We emphasize that by contrast to the flow of the cou- 
plings, for the quasiparticle weights Zj functions asso- 
ciated with the rcnormalizcd propagator, the whole flow 
must be computed with the fixed dressed Fermi surface, 
because the flow equations do depend on the Fermi sur- 
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FIG. 18: Flow of the norm J\f as s grows, corresponding to 
Fig. 1171 The inserted flow is a zoom on short times. 




1.5 1.55 1.6 1.65 



FIG. 20: Bare and dressed Fermi surface for A'^ = 32. 
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FIG. 21; Flow of the quasiparticle weights Zi = l/ipi, for 
A'^ = 32. We have represented these for J = 16 to 7 = 32, as 
the ones for smaller values of I can be found using the top- 
bottom {I ^ N — I because A'^ is even) symmetry. Note that 
we have indicated the value of t* = ln(Ao/A*) on the time 



face in the high-energy regime (see Eq. IjEMp ). As a 
consequence it is necessary to first compute the dressed 
Fermi surface, and then use it to compute the flow of the 
Zj^s. These flows are represented on Fig. 12 II We further- 
more show the variation of the Z/'s along the Fermi line 
at four different RG times on Fig.|211 We note that the 
dispersion in the .Z/'s is small, in qualitative agreement 
with the dynamical mean field theory results obtained by 
Biermann et al,^^ and this provides a consistency check 
of Bourbonnais' computation (see Sec. IIIIB3|l . Fig. 
shows that the evolution of the quasiparticle weights with 
typical energy scale exhibits some similarity with the re- 
sults of Kishine and Yonemitsu:''^ in the early stages of 
the flow, the quasiparticle weight is larger for ky = ±7r 
(/ = or iV), than for ky = 7r/2 {I = N/2), and this or- 
dering is reversed at later stages. However, we stress that 
the two models are different since Kishine and Yonemitsu 
have considered a 2D model with flat Fermi surface seg- 
ments, and it is not obvious that the end points of these 
segments should exhibit the same properties as the ex- 




FIG. 22: Here we show the transverse momentum (i.e. chain 
number I) dependence of the quasiparticle weights Zj, at four 
RG times t = 1, 2, 3 and 4 (see Fig. l21l for a time reference). 



FIG. 23; Flow of as a function of the "true" RG time 

t, for different values of A'^, and for = 0.1. These flows 

were obtained assuming that all coupling remain constant to 
their bare values (G"^ — —G^ — 0.3). Note that we chose 
/ = N/4: because is about the mean value of the Zj's. 

The inserted graph is a base 10 log-log representation of the 
value of the time derivative of — ln(Zjv/4) at time f = 10 as 
a function of A'^. The solid curve is the numerical fit, which 
shows a 1/A'^ behavior. 



tremal points ky = ±7r in our quasi ID model. The vari- 
ation of the quasiparticle weight along the Fermi surface 
has also been investigated by D. Zanchi^ for the 2D Hub- 
bard model, where he found a much stronger reduction 
of the Z factor in the vicinity of the Van-Hove singular- 
ities than for typical Fermi surface points. We believe 
this effect requires to take into account the variation of 
the Fermi velocity along the Fermi surface which we have 
not done here. The influence of these variations on the 
flow of couplings for a iV-leg Hubbard ladder has been re- 
cently studied^SiiS (for t± < t), with the conclusion that 
they play a dramatic role only below a cross-over scale 
which is extremely small as N becomes large. 

The second regime (15 < s < 60) is a transient between 
the ID high-energy flow, and the low-energy regime, 
where the shape of the Fermi surface is felt, and where 
the differentiation between the couplings takes place. In 
more physical terms, it corresponds to a Fermi liquid 
regime, located between a Luttinger liquid state at higher 
energies, and an ordered phase at lower energies. One 
might have expected that in this Fermi liquid regime, the 
Z/'s would remain constant, so that if this Fermi liquid 
regime was the final one, the quasiparticle residue would 
be finite. Instead of this, we see on Fig.|^ that the time 
derivative of the Z/'s decreases (in absolute value) but 
does not vanish. This is in fact a finite TV effect, as can 
be seen on the flows of the Z/'s one obtains (see Fig. l23|l . 
assuming that the couplings are constant, equal to their 
bare value, all along the flow. The inserted graph shows 
that the growth rate of the logarithm of the Z/'s behaves 
like 1/N in the Fermi liquid regime. This is indeed con- 
sistent with the flow equations (|E14p in the low-energy 
regime, where only the terms with a vanishing Aka{I, J) 
contribute, and whose number is of order N. This fac- 
tor N combined with the 1 /N'^ denominator explains the 
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numerical result. The relevance of these considerations 
with constant couplings is demonstrated by Figs. [THl and 
Indeed they show that the norm is almost constant 
for the interval t G [0, 5] in which we are interested (see 
Fig. EH). (Note that Figs. [H and [H were obtained for 
N ~ 8 not for = 32, but we have checked that on this 
interval t G [0, 5], the flows of the two norms are identi- 
cal, apart from a multiplicative factor of 8 = (32/8)'^/^, 
whose origin is the number N'^ of couplings for a given 
A'^, and which is irrelevant for our discussion). 

The final regime is one of a fixed direction, for which 
some normalized couplings are zero, and the others are 
gathered around specific values. It is in this final phase of 
the flow that the norm of the coupHngs explodes, as can 
be seen on Fig.^1 In order to be complete, we have also 
represented the link between the two RG times s and t 
on Fig.^l Notice that because of the definition of s (see 
the comment after Eq. (10) of our previous paper), and 
because the norm explodes in the end of the flow, the 
time t "saturates" , as a function of s, to a value which is 
roughly the critical temperature at which the final phase 
sets in. 

Let us now study more precisely the final fixed direc- 
tion, in the spirit of our previous paper. First of all, let 
us have a more precise look at the values of the couplings, 
on the fixed direction that is reached. These values are 
shown on Fig. [31 for the iV = 16 case. We did not choose 
= 8 as on Fig. El because we wanted to have more 
values (which was manageable here since we represent 
the whole set of values only once) . 

When briefly looking at these values, one can deduce 
that the couplings seem to be grouped into a few sets 
of similar values (with lots of couplings being equal to 
zero). Furthermore, forgetting about the zero value, it 
seems the three values (for charge or spin couplings) are 
not independent, but one is the sum of the other two. 
Finally, the values of the charge couplings seem to differ 
by a factor of three (and a minus sign) from the spin 
couplings. In fact, if we also look at Fig. 1171 we see that 
this will probably not be an exact statement for all values 
of N, but only in the limit of infinite N. 

It is then interesting to study what types of couplings 
take non-zero values. For the system we study, the nota- 
tion Gs{I , J) which is best adapted to superconductivity. 



can favorably be changed for F^''-^\l, J) — G'j^^j_g{I, J). 
5 is then the transferred transverse momentum, between 
the R particle that is destroyed, and the L particle that is 
created. In this notation, only F couplings with 5 = J —I 
or with 5 = N/2 are numerically found to have non-zero 
values. Notice that here, as N is even, N/2 is an inte- 
ger (we will discuss the odd N case a bit further). The 
couplings for which 5 = J — I and J — I ^ N/2 will be 
denoted as C'^'^^^ couplings, and correspond to the charge 



(respectively spin) couplings that are negative (respec- 
tively positive) on the fixed direction. They are the usual 
forward scattering couplings. The couplings that satisfy 
S = N/2 and J - I ^ N/2 will be denoted as V<^\ 
whereas the ones for which S = N/2 and J — I = N/2 
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FIG. 24: Values of the charge and spin couplings, for the fixed 
direction that is finally reached on Fig. 1171 but for = 16 
here. The inserted graph is a zoom on the values taken by 
the positive charge couplings, and shows that these can be 
grouped into two sets. 



will be denoted as A'^^^'^ couplings. Both have a trans- 
ferred transverse momentum S which is half the number 
of chains (i.e. tt if we use the usual momentum units). 
The vector linking a point of the Fermi surface on the R 
side, to the one on the L side, and N/2 chains further is 
a nesting vector, which explains why these couplings are 
present in the final low-energy fixed direction. 

We let the reader write down the RG equations sat- 
isfied by the F couplings, specialize these for the three 
types of couplings above, and deduce the equations sat- 
isfied for the final fixed direction, in the spirit of our 
previous paper. The resulting equations are: 
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FIG. 25: Flow of the charge and spin couplings, for the same 
values of the parameters as in Fig. [T7| but with iV = 7. 
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4^*^^ + 2(iV- 



= 4C''^ + 2(X'"2 - V^V") 
NV = {N - 2)(P'=2 + 3X>"^) + 



2(^'= - C'=)X>'= - 



(62) 



+ 2{A'' - C^)!?*^ 



4(^'^ + C'')V^ 



This set of coupled equations is nothing but Eq. (46) 
of our previous paper, with usual letters replaced by cal- 
ligraphic letters. The condition A"^'^ — C^'^ + D^'^ was 
satisfied, and ensured the SU(A^) symmetry of the in- 
teraction Hamiltonian. Here we will thus also be able 
to fulfiU the relation A''^''^ = C^""^ + V<^\ which was 
previously guessed when looking at the fixed direction 
obtained numerically. The values of the couplings can be 
found in Table II of our previous paper. It is clear that 
in the present situation, it is the so called (-I-, — ) fixed 
direction that is selected since in the infinite N limit, it is 
the one for which the charge couplings equal minus three 
times the spin couplings. 

The effective low-energy interaction Hamiltonian has 
the following schematic form (we drop the charge and 
spin structure): 



H 



cff 



9 I,k 
X 



cE[E4,/(fc+9)%./(fc) 

J.k' 

E [E4,7-Ar/2(^ + 9)cL„/(fc') 
J 

[E^L,/+W/2(^'"9)Cr,/(^) 



(63) 



X 



Let us describe the physics associated with such an ef- 
fective Hamiltonian. We will assume that N is large, so 
that we can neglect all finite N corrections. Thus, for 
example, only the T) terms (Peierls couplings) survive, as 
the forward couplings C arc a correction of order 1/iV. 
Furthermore, in the infinite N limit, the couplings take 
the values V = 3/A and 2?" = -1/4, so that V = -3X»". 
This relation implies that the interaction exists only in 
the triplet channel, and the effective Hamiltonian can be 
written as (g > 0): 
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.J-N/2,p 



{x)(Tp^p,lP^^j , (x) 



(64) 



E^L,/+A'/2,r(^)'^^.^"/'R,/,r'(2^) 



The only difference (except for changes of notation) 



between this effective Hamiltonian and the one we arrived 
at in Eq. (70) of our previous paper, is in the shift of the 
creation operators' chain number by an amount of N /2. 
We thus expect the physics to be essentially the same 
as we had discussed in our previous paper, apart from 
a different SDW's wave vector that will now be (2fc, tt) 
(remember k is the average Fermi momentum). We refer 
the reader to our previous paper for details. 



We thus have shown that after a high-energy ID regime 
where the Fermi surface's shape is not felt in the flow 
of the couplings, and after the crossing of the typical 
energy-scale given by the curvature of the Fermi surface, 
the system goes to a strong coupling phase of the SDW 
type, with the above effective low-energy Hamiltonian. 
The (2fc, tt) nesting vector naturally arises from the RG 
flow, and there is no need to artificially introduce it. 



Let us make a final remark, about the odd N case. The 
flow of the charge and spin couplings for A'^ = 7 is shown 
on Fig. 123 This figure is obviously different from Fig. 1171 
The reason for this is that there is no exact nesting vector 
anymore when N is odd. We have analysed which cou- 
plings are non-zero in the low-energy phase of Fig. 1251 
and these turn out to be BCS type couplings, indicating 
a superconducting low-energy phase. This is not in con- 
tradiction with what has been said before in the even N 
case. When N grows, the nesting is better and better in 
the odd N case, so that the RG flow will first be towards 
the same flxed direction as in the even N case. Then, 
there will be a shift from this fixed direction to another 
one, corresponding to superconductivity. But, this will 
take place at very low energies, and in regimes where the 
norm of the couplings has exploded. The conclusion is 
that the low-energy phase, in the thermodynamical limit, 
is always the one we have observed in the even N case. 
This discussion has been quite brief, but we refer the 
reader to our previous paper where we had analysed in 
detail how the observed shift from one fixed direction to 
another one, in a finite N situation, slows down as N 
increases and finally disappears in the infinite N limit. 
We have checked all this numerically, but unfortunately 
it requires quite a large value of N (more than 30) to be 
visible, so that we could not depict it in this paper. 
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Taking account of the Umklapps and limitations 
of the method 



In the above two simple examples we have studied, we 
encountered no real limitation of the computation scheme 
we have proposed. Of course, we had to use a non rigor- 
ous (but plausible) argument to define the scale at which 
we had to stop the flow of the couplings. We will see 
that there are some cases where it is not possible to use 
such a simple point of view. The well known main limi- 
tation of a RG approach is its perturbative nature. We 
cannot fully trust the RG flows when they go to strong 
couplings, even if the fixed direction that is reached in 
this regime gives an insight of what physics takes place. 
For the two computations of the Fermi surface we have 
given previously, it is clear that this was no limitation, 
because the deformation of the Fermi surface occurred 
in a weak coupling regime (remember Figs. 1171 and 1181 
where the norm diminishes between ^ s < 20, which 
is the time interval where the deformation of the Fermi 
surface takes place). 

The half-filled system is interesting, because it exhibits 
a variety of behaviors, depending on the strength of the 
bare couplings (compared to the value of the bare trans- 
verse hopping). We will discuss the strong, intermediate 
and weak coupling situations, which do not give the same 
low-energy physics. After this, we will consider the case 
of a nearly half-filled system. 



1. Half-filled system in strong coupling 

By strong coupling, we mean the initial couplings are 
large enough for the behavior of the system to remain 
purely one-dimensional. To make this statement more 
precise, let us study the flows for one of these strong 
coupling initial conditions. We will assume, as we al- 
ways did, that tj_/t^^ = 0.1. In this case, the Hubbard- 
condition G"^ = 0.2 = — and U = 0.4 is a strong 
coupling condition, for which the RG leads to a fully flat 
Fermi surface. In fact, as the couplings are large and 
grow quickly (because of the Umklapps), the Fermi sur- 
face flattens quickly, so that the decreasing cut-off never 
catches the scale of the Fermi surface, and there is no 
non-zero value of A*. The flow of the couplings is thus 
purely ID all along the flow and is well known, so that 
we do not show it. However, to be concrete, we show the 
evolution of the Fermi momenta on Fig. [211 About the 
norm, let us say that its value at the beginning of the 
flow is about 30, and at s = 140, it is about 1400, so that 
it is around 50 times bigger. This means the couplings 
have grown a lot. For example for the Umklapps, [/ ~ 19 
which is very big, so that the RG is not valid anymore. 
However, if we believe the RG is qualitatively valid, the 
flow of the Fermi momenta seems to show the existence 
of a conflned phase (the effective t± is zero). As the sys- 
tem goes to strong coupling, and remains purely ID, it 
would thus be natural to directly start from a system 
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FIG. 26: Flow of the Fermi momenta, for N = 16 chains, 
tx/i|| =0.1 and Hubbard initial condition = 0.2 = — 
and U = 0.4. 
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27: 


Flow of the three types of 


couplings. 


for iV = 



chains, — 0.1, and the initial condition — 0.03 

-G= and U = 0.06. 



of decoupled (no hopping) Luttinger liquids or ID Mott 
insulators. We shall simply direct the reader to some 
papers on this line of approach iii*^2i^ 



2. Half-filled system m weak coupling 

A weak coupling initial condition is one for which the 
Fermi surface does not get completely flat during the RG 
flow (i.e. a non-zero value of A* exists), and for which all 
couplings that are irrelevant (i.e. that do not exist at 
zero energy) do go to zero (after dividing by the norm) 
during the flow. As a consequence, the cross-over scale 
A* between the Luttinger and the Fermi liquid behav- 
ior is much larger than the typical scale for the onset 
of long range order. The system is therefore in a de- 
confined regime, in the sense that it allows for coherent 
transverse motion of electron-like excitations. An exam- 
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pie of this is obtained while using initial Hubbard cou- 
plings G"= = 0.03 = -G" and U = 0.06, and as usual 
t±/t\\ =0.1. It is not worth representing the deforma- 
tion of the Fermi surface in this case, for it is very small 
(for iV = 8, the effective is about 0.0995, so that 

the correction is of the order of half a percent). Let us 
however represent the flow of the couplings, on Fig. 1271 in 
the N = 8 case, and the flow of the quasiparticle weight 
(and of the norm of the couplings) on Fig. As in the 
incommensurate case, the time derivative of the quasi- 
particle weights becomes smaller (in absolute value) in 
an intermediate regime, and this effect is more and more 
visible as N gets bigger. This decrease can be understood 
from the flow of the norm, which shows a tendency to- 
wards a plateau behavior at intermediate scales, so that 
the arguments previously given in Sec. IIV Bl still apply. 
For the comparison with Fig. [57| let us simply say that 
for ^ t ^ 6 the link between s and t is approximately 
linear, and for t = 6, one has s = 14. 

Before studying the fixed direction, let us make a re- 
mark about the scale A*. Because of the presence of 
the Umklapps, not only the scales defined by all the 
Akp^a{I, J) play a role, but also the scales Afcp^(/, J). 
However, as the filling is one-half, the average Fermi mo- 



mentum is 7r/2, and one can check that in this case, the 
biggest AfcpQ,(/, J) is, as the biggest AfcF,a(-f, ^), equal 
to twice the difference between the biggest and the small- 
est Fermi momenta. A* will then be defined exactly as 
we did when the Umklapps were zero. 

The couplings on the fixed direction have very well 
defined values, with many being equal to zero. For 
the charge and spin couplings, the only non-zero cou- 
plings are the same as previously discussed in Sec. IIVBI 
namely of the type A, C and V. The condition A'^^^^ = 
Cc(s) j)c(s) gi^jjj seems valid, and one can furthermore 
check that C = so that A" = V. For the study 
of the Umklapps, it is also interesting to introduce V 
couplings, which are the equivalent of the F couplings: 
Vs{I,J) — Uj-i-s{I,J). An analysis of the non-zero 
Umklapps reveals that there are only three types of such 
couplings: U = Un,2{I,J ^ /), V = VV/2(/,J ^ /) 
and W — Un/2{I,I) — Vn/2{I,T)- Furthermore, these 
Umklapps also seem not to be independent, but linked 
by the W = 14 + V relation. We let the reader check 
that the fixed direction is found by solving the following 
set of coupled equations (which is just a generalization 
of Eq. 



' NA'' = (iV - l)iV''^ + 3P"2 +W + V^ - UV) 
NA"" = 4yt"^ + 2{N - 1) [V"^ + P^P'^ + \{pi^ - UV)] 
NC" = -(P'^^ -I- 3D''^ +UV-U^ ~ V^) 
iVC" = 4C"^ + 2 [V^^ - -f i(V2 - UV)] 
< NV = {N -2){V''^ + 3V^'^ +U'^ + V^ -UV) + 2{A'' -C'')V'' + 6{A'' ~C'')V'' + W{U + V) 

NV" = 2{N - 2) [X>"2 + V^V" + - UV)] + 2(^^ - C")X''= + 2{A'' - 0")^- + A{A^ + C")P'' + W{U - V) 
NU = 2{N - 2) [{V + 3P")W - 2P"V] + 2{V'= + P")>V + 2(^'= + 3^")^ - AA^V + 2(C'= - C'')U 
NV = 2{N - 2X2?'= - D")V + 2(2?'= - P")>V + 2{A'' - ^")V - 4C"W + 2(C'= + 3C'^)V 
NW = 2{N - 1) [ViU + V) + 3V%U - V)] + AA'^W 

. m 



It is easy to solve this set of equations (with the rela- 
tions between the couplings), order by order in N. One 
finds a few fixed directions, but the one of interest is the 
following (that we give to order 3, for the independent 
couplings): 
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We let the reader check that even the order reproduces 
quite accurately the values of the couplings of Fig. TT7\ 
(of course, up to an overall normalization factor). Here 
again, the effective Hamiltonian contains forward interac- 
tions that are 1/N corrections, and interactions for which 
the transferred momentum is the nesting vector. The 
2? couplings satisfy the relation + 32?** = (in the 
infinite N limit), so that the part of the effective low- 
energy Hamiltonian containing the 2? couplings is the 
same (apart from a numerical factor) as the one we pre- 
viously obtained (see Sec. IIV 51) . It is non-zero only in 
the triplet channel (we again consider the particle-hole 
parametrization of the couplings). Let us see what form 
the effective Umklapp Hamiltonian takes in this singlet 
and triplet parametrization. From Eq. H66(l . we see that 
we have a relation between U and V which reads V = 2U 
(this is valid up to 0{1/N^) terms). The corresponding 
interaction involving chains I and J on one side of the 
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FIG. 28: The top figure represents the evolution of .^jv/4 as a 
function of t for 7V=8 (squares), 16 (circles) and 32 (triangles). 
The bottom figure shows how the norm varies with time t. In 
fact, in order to allow the comparison, we have divided the 
norm for A'^ = 16 and 32 by a constant C so as to make all 
the norms equal at time t = 0. The squares thus represent 
the norm A/" for N = 8, the circles represent Af/C {N = 16) 
with C = (16/8)^/^ = 2V2 and the triangles represent JV/C 
{N = 32) with C = (32/8)^/^ = 8. 



Fermi surface, and 1 + 5 and J — J on the other side, may 
typically be written as (using Pauli's principle): 



2NL 



E E E^ii-.-'W-sir.p'W') 



I.J t.t' p,p' 



^ ('^RJ+W/2,T'^R,J-Af/2,p'^L,,/,p'CL,/,r' + ^- • (67) 



As Ir,T'Ip,p' — 2Ir,p'Ip,r' = —Ct.t' ' cTp.p' i it IS casy to 
rewrite the Umklapp interaction in the triplet channel 
only. If we define the generalized current 

(68) 

the total effective Hamiltonian takes the simple following 
form {g > 0): 
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(69) 



The low-energy physics can again be described by the 
fluctuations of the massless modes associated to the order 
parameter (which is ^{x) = {J^i^(x))). The difference 
with the non half-filled case studied in Sec. IIVBI is that 
the spin-density wave will be pinned to the lattice by the 
Umklapps. That this is indeed what happens can be seen 
by computing the effective action of the gapless modes, 
and one finds (we drop less relevant terms): 



N f 

Sesin) — — / dxdtdnnd^Ti. 



(70) 



with ^{x) = pn{x), p being a positive number found 
by solving mean-field equations, and n{x) is a real unit 



FIG. 29: Evolution of the smallest and largest charge 
(squares), spin (circles) and Umklapp (triangles) irrelevant 
couplings, for A'^ = 16, t±/t^^ — 0.1 and initial condition 
U = 0.02. 



vector, giving the direction of the staggered magnetiza- 
tion. This time, there is no gapless mode associated to 
the "phason" field (see the discussion around Eqs. (70) 
and (71) of our previous paper). This is physical, for a 
shift 9{x) 9{x) + 8 in the "phason" field corresponds 
roughly to a uniform translation of the spin-density wave 
condensate. As the physics of Eq. H70(l has already been 
discussed in our previous paper, we do not consider it 
further. 

Let us now study more precisely the fate of the irrele- 
vant couplings. We defined the initial coupling as a weak 
coupling if the normalized irrelevant couplings go to zero 
during the flow. In fact it is interesting to check that the 
flnal fixed direction that is reached is the same whether 
we run the complete fiow, or we run the fiow in which 
the irrelevant couplings are initially set to zero. We have 
done this in a system of TV = 16 chains, with t±/t^^ = 0.1 
and an initial Hubbard coupling U — 0.02. The evolution 
of the irrelevant couplings is shown on Fig. [23 In fact 
we have not represented all the couplings, because there 
are too many of them. We have decided to show only the 
smallest and largest charge, spin and Umklapp couplings. 
In order to make sure that the final fixed direction is the 
same as the one we would have obtained when initially 
setting irrelevant couplings to zero, we show the different 
values of the couplings on this fixed direction, in both 
cases, on Fig. [301 

A natural question that arises from this discussion is 
how small should the couplings be for being weak cou- 
plings according to the definition given above ? This 
question has already been answered by Lin, Balents and 
Fisher—. They have done so on a theoretical ground 
(and for a situation that is not the half-filled system, but 
this should not change anything), and found that the 
weak coupling condition reads U\n{N) ^ 1 (for large 
N). Thanks to our ability to take the irrelevant couplings 
into account, we have tried to check this numerically. To 
do so, we have determined the critical coupling Uc for 
which the irrelevant couplings do not fiow to zero any- 
more, for N =8, 12, 16, 24 and 32, and represented the 
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FIG. 30: Values of all the couplings on the final fixed di- 
rection, for — 16, t±/f|| = 0.1 and weak initial Hubbard 
coupling U = 0.02. On the left (respectively right), we rep- 
resented the values obtained when computing the whole fiow 
(respectively initially setting the irrelevant couplings to zero). 



values Uc In(A^) as functions of N. The result is shown 
on Fig. 1^ Because of the small values of N we have 
used, wc do not observe an horizontal line as could have 
been inferred from the U\n{N) <^ 1 criterion. But this 
criterion is in fact a sufficient condition (maybe not a 
necessary one) to observe a weak coupling behavior, since 
it implies that the effective Hamiltonian hardly changes 
during the high-energy part of the RG flow, for scales 
above A*. 



3. Half-filled system m intermediate coupling 

When the couplings are neither strong nor weak, that 
is intermediate, we suspect the system will behave more 
and more like a ID system, as the initial Hubbard cou- 
pling U grows. Before we check that this is the case, let 
us clarify this notion of intermediate coupling. We have 
just seen at the end of the previous section (jIV C 2|) . that 
the intermediate coupling should typically be character- 
ized by [/ C/c (remember Fig. In the case TV = 16 
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FIG. 31: Representation of the numerical values of Ucl'a{N) 
as functions of A'^, where Uc is the value above which the 
irrelevant couplings do not flow to zero. 




FIG. 32: Evolution of the (normalized) charge forward scat- 
tering couplings Go{I,J), when the strength of the initial 
Hubbard coupling U is increased {U = 0.12 for (a), (7 = 0.17 
for (b), U = 0.19 for (c) and U = 0.202 for (d)). The number 
of chains is A*' = 16, and tj_/t^^ — 0.1. The flow has been 
stopped at the RG time for which the biggest of the whole 
set of couplings is equal to 1. The charge forward scattering 
pictured here are the ones obtained at this time. In order 
to keep the flgures clear, we have not put any indication on 
the z axis. Let us simply say that the values of the couplings 
range from 3 • 10~^ to 1.3 • 10~^ in (a), and from 6.75 • 10" 



to 7.05 
(a). 



10 in (d), so that (d) is in reality much flatter than 




FIG. 33: Same as Fig. 1321 for the spin forward scattering 
couplings Go{I,J). Here, the values of the couplings range 
-2 • 10"* in (a), and from 7.95 ■ 10"* to 



from 



7.7- 10" 



S ■ 10"^ 
in (d). 



to 



and t_L/t|[ = 0.1 on which we shall focus, this means 
U ~ 0T5. We should also have U < 0.206, value above 
which the system is in the confined phase. If we ex- 
pect the system's behavior to change and become nearly 
one-dimensional for these typical values of U , this should 
mean that the effective hopping is of the same order of 
magnitude as the critical temperature. This will be dis- 
cussed in Sec. IIV C 41 when we study the phase diagram 
of the system. Let us simply say here that for the mini- 
mum (respectively maximum) value of the coupling U we 
will consider, namely U = 0.12 (respectively U = 0.202), 
the effective hopping — is about 8 times (respectively 
0.9 times) the critical temperature. These values confirm 
the previous expectation. 
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We have numerically studied the evolution of some spe- 
cial couplings, as U becomes bigger. It is not possible to 
consider the couplings on the final fixed direction. In- 
deed, the norm is huge even before it is reached. We had 
neglected this problem in the strong and weak coupling 
regimes. In the first case, we anyway knew that the RG 
is not valid anymore and should be replaced by a non- 
perturbative analysis. In the second case, we can make 
the norm as small as we want when reducing the initial 
coupling, because in this case, the RG fiow is scale in- 
variant (the RHS of the RG equations quickly become v 
independent, as v rapidly goes to very small values). In 
this last case we refer the reader to our previous paper— 
for more details about the implications of this scale in- 
variance. 

We thus have chosen to stop the RG flows at the time 
when the biggest of all couplings (the true couplings, not 
the normalized ones) reaches the value 1. In this regime, 
the RG should be valid (of course, the two-loops contri- 
butions are not negligible when the couplings approach 
1). This gives us the results shown on Figs. 1321 to 1511 
which were obtained for A'^ = 16 and initial couplings 
U =0.12, 0.17, 0.19 and 0.202. 

For the charge and spin couplings, we have represented 
forward scattering couplings Gq^*^' (/, J). In the weak 
couphng regime, we would have obtained Gq{I,J) ~ 
^j,i +N/2 (r emember the fixed direction we found in 
Sec. IIV C 2|l . Here we also obtain peaks around J = 
/ -I- N/2 values, but these peaks progressively disappear 
when U grows, as is expected because the system looks 
more and more one-dimensional. We have chosen to 
represent Us{I,I) couplings, in the case of the Umk- 
lapps. The reason for this choice is that the biggest 
of all Umklapps is (numerically) found in this subset of 
couplings. Again, the weak coupling would give a peak 
Us{I,I) 5s^is!/2, which is smeared in the case of inter- 
mediate couplings, and disappears in the strong coupling 
limit. 




FIG. 34: Same as Fig.^l for the Umklapp couplings Us{I, I). 
Here, the values of the couplings range from 8 ■ 10~^ to 2.4 • 
10"^ in (a), and from 1.3954 • 10~^ to 1.3961 • 10"^ in (d). 
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FIG. 35: Phase diagram of the system, computed for N = 8 
drains. See text for a detailed description of this diagram. 



4- Phase diagram 

As a conclusion of this investigation, we shall summa- 
rize the numerical results we obtained on a single figure, 
which is the phase diagram of the system. It is depicted 
in Fig. 123 

The solid curves represent as a function of t±/t^^ , 

for [7=0.05, 0.06 (indicated by the arrows), 0.08, 0.1, 0.15 
and 0.2 (notice that both x and y scales are expressed in 
base 10 logarithms). For a given U, is zero in the 
confined phase, for t± smaller than a critical value 
(this explains the vertical lines), and it takes non-zero 
values as soon as t± > i'^. When t± is much larger than 
t']_ , —t±, so that all solid curves asymptotically go to 
the first bisector. The dashed curves give the value of the 
scale at which the couplings diverge, which is the critical 
temperature Tc. They are horizontal when t± < i^, since 
in our approach the RG flows are purely one-dimensional 
in this regime. The upper (respectively lower) dotted 
curve is a straight line of slope 1 (numerically found 
to be 1.004), going through the points of coordinates 
(tliU) , Tc (U, tl {U)) (respectively {tl (U) , tf (tl (C/)+)) , 
as the one represented by a diamond (respectively circle) 
in the inserted figure. This inserted figure is a zoom of the 
interesting region where both scales meet, for U = 0.05. 
We have indicated the different phases (Luttingcr Liquid, 
Fermi Liquid, Mott Insulator and Spin-Density Wave). 
The dash-dotted curve is the first bisector, that we did 
not represent in the global figure to keep it readable. 

Let us now study the physical implications of the value 
1 taken by the slopes of the two dotted lines. The up- 
per one tells us that the critical value of t± is propor- 
tional to the value of the charge gap A of the Mott 
insulating phase, and we numerically find ~ 1.14 A. 
The lower curve gives the following relation between the 
value of the effective hopping for t± = ^1^, that we 
will denote by f^^* , and the critical value of the bare 
hopping: f^^* 0.68i'i^. Let us remark that this also 
implies: i]^* ~ 0.78A. These relations show that the 
confinement-deconfinement transition takes place when 
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the bare and the effective hopping become of the order 
of the charge gap. These results have natural interpre- 
tations. Imagine the system is in the Mott insulating 
phase, with zero effective hopping. It is clear that we 
can only compare the bare hopping to the charge gap. If 
on the contrary the system is in the deconfined regime, 
the low-energy physics is dominated by two scales: the 
critical temperature and the effective hopping. When de- 
creasing the bare hopping, we expect a phase transition 
to occur when these two low-energy scales become of the 
same order of magnitude. These results are in quantita- 
tive agreement with those obtained by Tsuchiizu et al^ 
for a two-chain model, see for instance the insert in their 
Fig. 1. For values of t± which are not too small, they 
indeed find a proportionality between A and t±_ with a 
slope compatible with our results. For smaller values of 
t±, they obtain a sizeable deviation away from a linear 
behavior. But this difference with our conclusions comes 
from their choice of a fixed value for the intra-chain for- 
ward scattering, while they let the Umklapp scattering go 
to zero. In this regime, they observed a significant renor- 
malization of the hopping amplitude, so that the transi- 
tion is finally given by the balance between the charge 
gap and the renormalized hopping. 



5. Nearly half-filled system 

Up to now, we have not encountered any real difficulty 
when choosing the A* scale. Of course, the choice was 
purely pragmatic, as we simply chose the biggest of the 
scales K appearing in the RHS of the RG flows. There 
was no problem in the half-filled case, because the biggest 
scale was the same for all the nine sorts of K (remember 
the definitions in Eqs. (|E4|) to IIE12|l ). Let us consider 
the non half-filled case, for which the mean Fermi mo- 
mentum is k = tt/2 + 6k. The filling does not change 
the biggest K^^ and K^^, which is still the difference 
between the biggest and the smallest Fermi momenta, 
Afcj?^'^ = fc^^'^ - kf'\ We let the reader check that the 
biggest if^^ is now 2\Sk\ + Ak^"""" and the biggest K'^'^' 
is \6k\ + Akp'^^. Those last two scales are obviously al- 
ways bigger than the first one. 

What are the consequences of the existence of these 
three scales? At the formal level there is no real conse- 
quence. Indeed, even in the half-filled case, there were 
a lot of different scales, given by all the possible K's, so 
that introducing more scales does not make much change. 
But, practically, we have used the pragmatic point of 
view that we should stop the cut-off scaling flow when 
the biggest scale is reached. This relied on the hypothe- 
sis that the flow of the Fermi surface nearly stops at that 
scale. It is clear that if the filling is not too far from 
one-half, the three scales are not qualitatively different, 
so that we can afford stopping the flow when the biggest 
scale is reached. 

When the filling is quite far from one-half, things be- 
come more involved. Of course, we could simply for- 



get about the Umklapps and perform the analysis of 
Sec. II V Bl devoted to the non half-filled system. But we 
expect that in an intermediate to strong coupling regime, 
the Umklapps could play a non-negligible role in the high- 
energy part of the flow, where there are not yet irrelevant. 
We would thus like to be able to take them into account. 
Intuitively, we expect that the Fermi surface deforma- 
tion is caused by both U and G couplings for a cut-off 
A > Afcp'*'' -I- {2)\Sk\ and by the G couphngs only for 
Afcp'"' < A < A/cp'"' + {2)\5k\. It is however not possi- 
ble to implement this idea in a simple manner. Indeed, 
once the scale Akp'^^ + (2)|^fc| is reached, we could drop 
the Umklapp contribution to the flow of the Fermi sur- 
face, but the flow of the Umklapp couplings will not stop 
and will depend on the shape of the fixed dressed Fermi 
surface. This flow of the Umklapps will affect the flow of 
the G couplings and as these latter still deform the Fermi 
surface, we see that the flow of the Umklapps indirectly 
affects the flow of the Fermi surface. As the shape of the 
dressed Fermi surface comes into play before the flow of 
the running Fermi surface stops, there is here no simple 
way to circumvent the inversion problem we have dis- 
cussed in the introduction. 

As a consequence, in what follows, we will restrict our- 
selves to situations where our pragmatic scheme works, 
i.e. to the nearly half-filled system. We have considered a 
fiUing slightly less than 1/2, setting the chemical poten- 
tial to -0.01. For this value, the difference between the 
initial values of Ak^""" and Afc™^^-^- (2)|(5fc| is about 10%, 
which is reasonable. As in the half-filled case, we have 
observed different regimes, when changing the strength 
of the initial coupling namely weak, intermediate and 
strong coupling regimes. 

In the strong coupling regime, as before, the system 
behaves as a purely ID system, and the Fermi surface gets 
completely flat. The intermediate coupling regime is the 
same as the one observed in the half-filled case, where 
the Umklapps are not suppressed, because the scale at 
which the phase transition takes place is bigger than the 
scale that measures the distance from half-filling, namely 
{2)dk. In the weak coupling regime, the Umklapps are 
irrelevant, and all vanish (for the normalized couplings) 
at low energies. The final fixed direction is simply the 
one we previously found in Sec. II V Bl where we had set 
the Umklapps to zero at the beginning of the flow. The 
flows of the couplings are represented on Fig. 1361 and 
were obtained for N = 8 and initial coupling U = 0.1. 
We also have represented the evolution of the norm of the 
couplings, because its behavior changes drastically at the 
precise time the Umklapps become irrelevant (s ~ 24). 
The plateau observed in the norm's evolution just after 
this time reveals the existence of an intermediate Fermi 
liquid phase. 

These flows show how our RG, taking account of the 
different scales of the system, is able to get rid of the 
irrelevant couplings, such as the Umklapps, and leads to 
the correct flnal fixed direction. This is not of a purely 
academic interest. Indeed, if one takes the same value of 
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FIG. 36: Flow of the three types of couphngs, for N = 8 
chains, tx/ty ~ 0.1, /i = —0.01, and the initial condition 
= 0.05 = and U = 0.1. We also have represented the 
flow of the norm, for ^ s ^ 60. 



the initial couplings {G^ — 0.05 = —G^), but set U to 
zero from the beginning (with a flow very much hke the 
one in Fig. I17|l . the scale at which the phase transition 
occurs is found to be 4 ■ 10^^ times smaller as the one 
found when the Umklapps are incorporated and vanish 
along the flow. This can in part be explained because ne- 
glecting the Umklapps from the beginning of the flow is a 
very crude approximation. We could also take the Umk- 
lapps into account in the ID part of the flow, and then 
take them to zero. However, a look at Fig.l36lshows that 
the Umklapps do not vanish very fast, and their influence 
has thus no reason to be small. Indeed, we performed this 
comparison, and found that this prediction is correct. In 
fact both methods neglecting the Umklapps at one time 
or another, give approximately the same critical temper- 
ature, because the Hubbard coupling U is small so that 
the couplings do not change much in the ID region (we 
found 6"= = 0.074 and ^ -0.038 in the end of the ID 
flow). 



V. CONCLUSION 

We have attempted to develop a simple physical pic- 
ture to understand the forces which drive the defor- 
mations of the Fermi surface of an interacting electron 
system. Using considerations from second-order time- 
independent perturbation theory, we showed that the 



shape of the dressed Fermi surface controls the quantum 
zero-point motion correction to the ground-state energy. 
We demonstrated that a given coupling tends to deform 
the Fermi surface so as to have all its four momenta 
precisely on the Fermi surface, because this allows for 
smaller energy denominators (in the second order contri- 
bution) and thus decreases the total energy. As a con- 
sequence, we find that the Fermi surface is deformed by 
irrelevant couplings, which are characterized by the im- 
possibility of choosing all four momenta on the Fermi 
surface. Because of this, the Fermi surface deformation 
induced by these couplings can only occur in the high- 
energy regime where the kinematical constraints associ- 
ated to the Fermi surface do not play much role. Once 
the energy is low enough and the warping of the Fermi 
surface is felt, these irrelevant couplings do not flow any 
more and no longer contribute to the Fermi surface de- 
formation. 

For the quasi ID materials at half-filling, the Umklapp 
couplings belong to this category. As they undergo a 
strong renormalization in the high-energy regime, they 
have a much more drastic effect than the charge or spin 
couplings in the doped system. Our numerical simu- 
lation provide a description of the cross-over from the 
confined regime to the Fermi liquid, which is in over- 
all good agreement with previous works AiSiSLiSiS We 
have presented a detailed analysis of the evolution of 
the quasiparticle weights as a function of the typical en- 
ergy scale. This confirms the existence of a Fermi liq- 
uid regime at intermediate energies, for the deconfined 
systems. It would be very interesting to compute the 
longitudinal and transverse optical conductivities, which 
could be done by adapting some existing methods^ to 
quasi ID systems. Another problem is to investigate the 
nature of spin correlations in the confined regime. This 
can not be achieved within the present formalism since 
the couplings diverge at a scale associated to the charge 
gap which is much larger than the Neel temperature in 
the limit of small transverse hopping. This problem has 
been addressed by Kishine and Yonemitsu^^ who used 
RG equations to two-loop order for the couplings. But it 
is not clear that the two-loop corrections provide a reli- 
able description since the couplings do not remain small. 

This limitation of our method is certainly connected 
to the fact that we are using a physical picture in which 
the fermion fields remain the elementary objects. This 
is valid at sufficiently high energies and thus generi- 
cally adapted to the study of the Fermi surface defor- 
mations. However in the confined regime, elementary 
excitations are likely to be very different from the Fermi 
liquid-like quasiparticles, but rather some soliton-like ob- 
jects. In this case, it seems a deeper understanding of 
the corresponding phases should be obtained by expand- 
ing around the exact solution for a system of uncoupled 
chains4i*^2*^ Nevertheless this raises the important is- 
sue of the validity of an adiabatic principle for gener- 
ating the ground-state at finite transverse hopping from 
the ground-state of uncoupled chains. At half-flUing this 
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adiabatical principle certainly holds in the confined phase 
where the charge gap is finite. But its validity beyond 
the critical value of the transverse hopping is question- 
able. This may explain the qualitative discrepancy be- 
tween the perturbation approach of Essler and Tsvelik;^ 
which leads to a disconnected Fermi surface with electron 
and hole pockets, and the dynamical mean-field theory of 
Biermann et al, who have obtained a conventional Fermi 
surface (see Fig. 5 of the latter work^). Away from half- 
filling, this notion of adiabatic continuity is even less ob- 
vious to prove since the energy gap vanishes for uncou- 
pled Luttinger liquids. However we believe the use of a 
skeleton expansion by Arrigoni'^^ likely provides a way to 
circumvent this potential difficulty. 
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APPENDIX A: EQUIVALENCE BETWEEN 
MINIMIZING THE TOTAL ENERGY AND 
STABILITY CRITERIA FOR THE DRESSED 
SINGLE PARTICLE PROPAGATOR 



1. Two-chain system at first order 

We wish to check that the self-energy computed by 
using a free excited state, with a trial Fermi sea which 
should be determined only at the end of the calculation, 
gives the correct results of Sec. lIIB'2l The sole difference 
with the calculation of Appendix IB II below is that we 
shall consider here the following free propagators: 



turbation theory) that the self-energies read: 



/i(0) + v^°\k - 4j) + i^sgn(fc - kpj) 



/i(0) - 4°^ (fc + 4"/) - i??sgn(fc + fcp 



, (Al) 



(A2) 



These formulas differ with (|B1I) and ljB2p below only 
in the imaginary parts. We let the reader check (some 
more details about one-loop self-energy calculations can 
be found in Appendix IB II devoted to the standard per- 



S«(fc,c.) = A 



Ao + {kF,o ~ 4'!o) 



27r 



Ao - (fcp.o - 4'!o) 



-C 



27r 



Ao - (fcF,0 - 4"o) 



27r 



Ao + (fcF.o - 4"o) 



-C 



27r 



(A3) 



(A4) 



It is sufficient to compare these results with (|12f) and (|13f) 
to understand that we will exactly find H14|l and H17|l . 
The energy minimization method or this self-consistent 
computation of the self-energy carried formally to first 
order in interaction do generate the same higher order 
terms. These appear since corresponding contributions 
are sensitive to the shape of the trial Fermi surface. As 
shown in Appendix ^ this feedback effect is missing in 
the standard perturbative approach. 



2. Formal proof for a finite size system 

Let us choose a trial Fermi surface, with the corre- 
sponding occupation numbers n{k) G {0,1} and $(fc) = 
1 — 2n{k). This generates a free particle state (Slater 
determinant) which may be used as a starting point 
for perturbative expansions of the total energy E({n}) 
and the single particle propagator G<j,{k,uj). Suppose 
we add one particle to the system, so that the total 
momentum is increased by kg. This is achieved by us- 
ing n'{k) = n{k) + 5k,ko instead of n(fc), assuming that 
^(^o) = 0. This induces a change in the total energy of 
the system AE — E{{n'}) — E{{n}). For a finite size 
system, it is easy to connect this energy shift to the sin- 
gle particle propagator G'$(fco,ti^). Indeed, we have the 
well-known spectral decomposition: 



G^{ko,uj) = ^ 



{N+l,ko,a\cl\N) 



uj - E{N + 1, /cq, a) + E{N) + iry 



\{N-l,-ko,P\c^^N)\ 



' + E{N -l,-ka,(3)-E{N)~iri 



(A5) 



Here, | A^) is the eigenstate with N particles obtained per- 
turbatively from the free-particle state with distribution 
n(fc), and the kets \M,k,a) denote eigenstates with M 
particles and a total momentum k with respect to the 
total momentum of state |A^). The total energy of these 
states is of course E{M, k, a). So the energy shift AE is 
one of the poles of Gg,{ko, uj) seen as a function of oj. De- 
noting the typical interaction strength by V, it is sensible 
to assume that for a finite-size system, energy differences 
like E{N + 1, ko, a) - E{N) (and AE in particular) can 
be expanded as power series in V. In the non-interacting 
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case, Ai? = £0(^0)1 and {kQ,uj)~^ = lu ~ £0(^0) + ^V- 
Therefore, AE is the pole of G'$(fco,a;) (as a function of 
u)) which goes smoothly towards £0(^0) as V goes to zero. 
Writing G$(fco, ti-")"^ = w — £o(A:o) — S$(A:o, w), and given 
the fact that S$(A:o7 '-^) has a well-defined power series ex- 
pansion in V which vanishes as V goes to zero, we may 
conclude that AE may be obtained as a formal power 
series in V from the solution of the following equation 
for uj: 



uj — £0(^0) ~ 5RS$(/co7 ^) ~ 0. 



(A6) 



Indeed, let us denote by w(fco) the solution of this equa- 
tion which goes to £0(^0) as V goes to zero. Then, we 
have: 



AE = uj{k()). 



(A7) 



Now, if the trial state | A^) obtained from n{k) minimizes 
the total energy £'({n}), it means that removing a parti- 
cle at ki on the Fermi surface associated to the distribu- 
tion n{k) and adding another particle at k2 also on the 
Fermi surface does not change the total energy (up to 
corrections which are negligible for very large systems). 
This yields w(fci) — ti'(fc2), implying that quasiparticle 
energies are constant (equal the the dressed chemical po- 
tential ^) on the Fermi surface associated to n{k). This 
is exactly condition i) (see Eq. |0J|) for the dressed propa- 
gator discussed in Sec. Ill Al Assuming condition i) holds, 
condition ii) on the imaginary part follows from standard 
phase-space arguments and analyticity considerations de- 
veloped already long ago by Luttinger— or Langer— . 
The main idea is to use an expression for the self-energy 
in terms of skeleton graphs. Condition i) suggests that 
the full one-particle spectral function (which determines 
completely the internal lines of these graphs) is qualita- 
tively similar to the one of a Fermi liquid with the Fermi 
surface obtained from n{k). 



3. Extension to an infinite system 

Applying the previous argument to infinite systems re- 
quires some care. In fact, we have to prove that the co- 
efficients of perturbative series in powers of V for AE 
or aj(fco) have a well-defined infinite volume limit. Our 
experience with other systems including models for an 
unstable state coupled to a continuum suggests that such 
a limit does not exist in general. However, perturbation 
theory in powers of the interaction strength for fermion 
systems with local two-body interactions is likely to be a 
favorable case for which this limit may be safely taken. 
For instance, the perturbative expansion of the ground- 
state energy involves connected Feynman graphs with no 
external lines, which contributions are easily shown to 
be proportional to the volume. Similarly, standard tech- 
niques based on the Luttinger-Ward energy functional 
(see for instance the text by Nozieresji^ pages 222 to 



229) show that: 
AE ^ E{{n'}) - E{{n}) 

= £o{ko) - 



TT^log 



w - soiko) - i?? 



- £o(fco) + i?7 
- £o(fco) - S$(fco,c^) - 



(A8) 



177 



^ - eo(fco) - S*(fco,t^) - ir] 



up to terms which vanish in the thermodynamical limit. 
If fco is close enough to the dressed Fermi surface so 
that the inverse life time of the corresponding "quasi- 
particle" is small compared to 77, it is easy to show that 
AE — uj{ko). This shows that the series expansion 
of (jj(fco) has a well-defined infinite volume limit. This 
fact is a priori non trivial since any perturbative algo- 
rithm for uj{ko) involves partial derivatives at any order 
for J^^{ko,uj) with respect to to, taken at w = £o(fco). 
Although I]$(fco,w) has a good thermodynamical limit, 
some difficulties arise while considering derivatives with 
respect to to. Indeed, their expressions for finite size sys- 
tems involve sums of rational functions of lo with multiple 
poles, and these are not easily converted into converg- 
ing integrals in the infinite volume limit. But the above 
connection between AE (which has a thermodynamical 
limit) and uj{ko) shows that all the wild terms which are 
expected to appear in a perturbative expression of ^(fco), 
eventually cancel. 



APPENDIX B: DIFFICULTIES WITH THE REAL 
PART OF E IN THE TRADITIONAL 
PERTURBATION SCHEME 

1. Two-chain system at first order 

Let us begin by the calculation of the self-energy in the 
usual case where one starts from the free ground-state. 
The free propagators are simply given by: 

.(0) 



.(0) 



(Bl) 



(0) 



v^°Hk + k^"j) - ifjsgni k + k'^'j) 



(0) 



-. (B2) 



We will restrict ourselves to the study of the right propa- 
gators, because the left ones can be analyzed in an analo- 
gous way. The first order correction to the right propaga- 
tors is given by the tadpole graph represented on Fig. 1101 
where the solid (respectively dashed) lines represent right 
(respectively left) propagators, and where the black dot 
denotes one of the couplings. The self-energy for the 
(R,0) fermions is given by two terms: either the interac- 
tion is j4, in which case the left propagator in the loop is 
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on branch 0, or the interaction is C, and the left propaga- 
tor is on branch tt. It is a simple matter to evaluate the 
tadpole, and to show that in the thermodynamic limit 
the self-energy given by the A interaction reads: 



dq 



(B3) 



free propagators corresponding to the dressed Fermi sur- 
face. A similar conclusion also holds for the chemical 
potential shift, as a comparison between Eqs. (|17|l and 
(|B9|) readily shows. 



Formal calculation to second order 



where «l°q('7) is the particle distribution on branch (L,0), 
i.e. it is 1 if g ^ — fcp.o and otherwise. Of course, ex- 
actly as in the energy minimization scheme described in 
Sec. IIIBI we get infinite results because our linearized 
dispersion relations have been extended to include in- 
finitely many states. We will thus here too regularize 
these divergences by putting an ultra-violet cut-off Aq 
on the momenta, around the four free Fermi momenta 
(remember Fig. 0] for one band). It is then easy to show 



that E 



(1) 

R,0;A 



ylAo/(27r). We let the reader check that 
the final results for the self-energies of right fermions are: 



(B4) 
(B5) 



The renormalized chemical potential fj, and Fermi mo- 
menta fcF,o(ir) can now be deduced from the condition 
that the inverse propagators vanishes for uj = fj, and 
k = fcp.o or fc = /cf.ttj and from the conservation of the 
number of particles. This last condition is nothing but 
the Luttinger theorem. We thus have to solve for the 
following system of three equations for three unknown 
quantities: 



,(0) 



,(0) 



■ V 



f (fcp.O 
(0)/ 



,(0) 



{A + C)^^Q, (B6) 



A*^-+<^(fcF,.-A:n) -{B + C)^ = 0, (B7) 

ZTT 



«F,0 



(0) 
F.O 



k 



(0) 
F,7r- 



(B8) 



The chemical potential is found by summing the first two 
equations and making use of the third one. Then one 
gets the difference between interacting and free Fermi 
momenta at one loop: 



,(1) 



An 
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^F.O ~ ^F,0 ^ ^ ^)" 



Ao 



(0) 



(B9) 
(BIO) 



These are the results given by the standard perturba- 
tion theory. This last result is to be compared with (|14|l . 
To first order in the couplings, both results are equal. 
But (|14l) contains next order contributions that are not 
present in (|B10|) . This happens although both compu- 
tations assume the same physics, namely the validity of 
the Hartree approximation. As has been shown in Ap- 
pendix lA II consistency between both viewpoints is re- 
covered only if the electron self-energy is computed with 



Let us consider a system of interacting spinless 
Fermions in c? = 3 dimensions. In fact, the actual value 
of d does not have much infiuence in the following discus- 
sion, the main point is that d ^ 2 so the Fermi surface 
is in general a smooth manifold of codimension one in 
fc-space. We take the following Hamiltonian: 



H 



V r d^k 
2 



e{k)c1{k)c{k) 
d^k' 



(Bll) 



d^q 



/(fe,fe',g) 



(27r)3 J (27r)3 J (27r)3 ' 

xct(fc + 9)ct(fc' - q)c{k')c{k). 



We assume the Fermi surface for V = is connected and 
that each half-line starting from the origin in fc-space in- 
tersects it only once. For any unit vector tt, wc thus 
define a positive number fcF,o(^*) such that kp_o{u)u be- 
longs to the Fermi surface. The Fermi sea is then the 
set of k points such that fc = ku with u unit vec- 
tor and ^ k ^ kpfllu). The total particle number 
is assumed to be fixed, independently of the coupling 
strength V. We now consider eigenstates obtained by 
adiabatic switching of the interaction V on free particle 
states with a deformed Fermi surface u kF{u). De- 
noting kp(u) — fcF.o(it) = Skpiu), the constraint on the 
total particle number reads: 



d^u 



kl f^{u)SkF{u) 



(B12) 



1 



+kF^oiu)Skp{u) + -Sk^{u 



= 0, 



where d^it is the usual area element on the unit sphere, 
for instance d^u — sm0d9d(j) in spherical coordinates. 
We now wish to choose dk-p{u) in order to minimize the 
total energy of the corresponding eigenstate, while keep- 
ing a constant particle number. To second order in V, 
this total energy is given by: 



mkF}) 

V 
+ 2 



^-(fc)e(fc) 



(B13) 



d^fc 
(2^ 
d^fc 



d^fe^ 

d3fc' 



n{k)n{k')g{k,k',0) 

(27r)3 J (2^)3 J (2^^^^'^''^)' 
n(fc)n(fc')(l - n{k + q)) (l - n{k' - q)) 



e{k) + e{k') - e{k + q) - e{k' - q) 
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Here, n{ku) = 1 if fc is smaller than k-p{u) and 
n{ku) = otherwise. We have also defined g{k, k', q) = 
f{k, fc', q) — f{k, k', k' ~k~ q). After some algebra, we 
find that 5kF\u) = V5kY^i{u) + V'^5kY,2{u) + 0{V^), and 
^ = ^iQ + Viii-^V'^^i2 + 0{V^), where 5k'p^i{u), 5fcF,2(M), 
/ii and /i2 are given by the following coupled linear equa- 
tions: 

VF{u)5kF^i{u) + E(i) (fep^oH) - Ml = 0, (B14) 



y" d^Mfcl^oH'^^F^H = 0, (B15) 



fc^_o(«)^feF,2(M) +fcF,o(M)<5fc|,l(M) =0. (B17) 



In these expressions, vy{u) and v'^iu) denote the first 
and second derivatives of the function x t-> e((fcF,o(''^) + 
taken at x = 0. T,''^\k) and S(^'(A;,cj) are the self- 
energies computed to first and second order in V ^ using 
the standard algorithm: 



+u ■ VfeE(i)(fcF,o(M)) + ^v'p{u)Skl^^iu) (B16) 
+AE(i)(fcF,o(M)) - Ai2 = 0, and 



d^fc^ 

(2^ 



i^"^k')g{k,k',0), and (B18) 



«^(2)/, A I f d^k f d^k' f d^q , 



(B19) 



(27r)3 ./ (27r)3 J (27r)3 

n(")(fc')(l - n^^Kk + q))(l - n(")(fc' - q)) + (l - (fc'))f^^°Hfc + qW°Kk' ~ q) 



uj + e(fc') - £(fc + q) - £(fe' - q) 
I 



The most interesting quantity in these formulae is 
AE(^^(fe). It is the change in the first-order (with re- 
spect to V) self-energy due to the fact that the Fermi 
surface has changed by an amount (5fcF,i- More precisely, 
we have: 



d^u' 



g{k, fcF.o(w'), 0)fcF,o('"')'5fcF,i(w')- (B20) 



It turns out this term is not recovered in the naive per- 
turbation algorithm. This latter procedure is based on 
solving for kp(u) in the equations: 

£(A;f(m)) + 5ftE(fcF(M), = /i, (B21) 

where E(fc,a;) is computed with the free propagators as- 
sociated to the non-interacting Fermi surface. As before, 
/i is chosen to keep a constant total particle number. 
This second approach yields the same set of equations 
as before, except that the term AE^^^ (feF,o(''^)) is miss- 
ing in the first equation for Skp^2{u). This shows that the 
naive algorithm is not able to keep track of the first-order 
Fermi surface deformation while evaluating the Hartree- 
Fock corrections to second order. Intuitively, these effects 
are expected to be associated to the four second order 
graphs for S shown on Fig. 1371 However, these graphs are 
anomalous according to Kohn and Luttinger, and their 



contribution vanishes in the T = perturbation theory 
scheme. We believe this illustrates the crucial problem 
with naive perturbation theory. As we go to higher or- 
ders in V, lower order graphs for E are modified by the 
changes already induced on the occupation numbers of 
single particle states. One would expect to capture these 
changes thanks to self-energy insertions in the internal 




FIG. 37: The four anomalous Kohn-Luttinger graphs con- 
tributing to the self-energy at two loops. 



32 



lines of the lower-order graphs for E. But the example 
of anomalous graphs shows this does not work so well in 
general. 

As discussed in section III CI the natural cure for this 
problem is to fix the dressed Fermi surface, thanks to 
counterterms which gradually modify the single particle 
dispersion of the original Hamiltonian, as V is increased. 
In this approach, we choose the dressed Fermi surface 
u 1-^ kF{u) and the dressed dispersion relation e(fe). We 
therefore compute the self-energy l^{k,Lu) with respect 
to this dressed Fermi surface. Denoting by = e{kp{u)) 
the dressed chemical potential, the counterterms "EcTiu) 
are defined by: 

E{kF{u),fi) + ^CT{u) = 0. (B22) 

In this third approach, the "bare" Fermi surface u i-^ 
kpfliu) becomes a function of V and {fcp}. It is obtained 
from: 

e{kF,o{u)) + Sct(m) = MO, (B23) 

where, as always, /ip is chosen in order to conserve the 
total particle number. We therefore have to solve: 

e{kF{u)) - e{kFfl{u)) + I]{kF{u),fi) ^ fi - fiQ. (B24) 



Since S is computed with free propagators whose sin- 
gularities lie on the dressed Fermi surface, it is easy to 
check that this yields the same expressions for 6kF,i and 
(5fcF,2 as the energy minimizing procedure, in complete 
agreement with the general conclusions of Appendix 1X1 



APPENDIX C: DIFFICULTIES WITH THE 

IMAGINARY PART OF E IN THE 
TRADITIONAL PERTURBATION SCHEME 



Let us consider the traditional perturbation scheme 
around the unperturbed Fermi surface. With the same 
notation as before, this corresponds to the choice of free 
propagator: G(°)(fc, cj)^-^ ^ uj — eo{k) + irisgn{eo{k) — 
fio). In the discussion, we shall use the spectral densities 
Pp^h{k,u)) for excited states involving p particles and h 
holes, with a total momentum k and a total energy ld. 
We have: 



Pp,H{kM = f{J (fp^(so(fe.)-Mo)n/ (^^(M0-£0(fc;)) 



dfc; 

p h p h 

i—1 j — 1 i—1 j — 1 



(Cl) 



Let us consider the simple second order diagram for 
E(fc,ti;) shown on Fig. ^] (the sunrise diagram). It is 
simple to check that the imaginary part of this diagram 
is proportional to p2,i{k,Ljj) — pi2(— fc,— w). As is well 
known since Landau, this quantity vanishes at the bare 
chemical potential /io, and behaves as (w — pof' in mag- 
nitude for Lo close to /ip. The effect of Fermi surface de- 
formation on arises via the replacement of bare 
propagators by sequences of these propagators separated 
by lower order self-energy insertions. The main point we 
wish to emphasize here is that there is no simple way 
to predict the influence of these insertions on the fre- 
quency dependence of For some graphs, and 
some patterns of insertions, the resulting 3E(fc,Ci;) will 
continue to vanish at w = /xq, whereas some other com- 
binations will produce a finite contribution to '^'il(k,uj) 
at w = /LtQ. Therefore, the traditional scheme does not 
allow for a good control of the analytical structure of the 
self-energy. 



Let us show this on a typical example. The simplest 
interesting situation is obtained for the sunrise graph, 
for self-energy insertions which are assumed not to de- 
pend on k nor on lo, since it is then easy to perform the 
frequency integrals. In the more general case of an ar- 
bitrary frequency dependence for the insertions, the nat- 
ural procedure would be to Taylor-expand them in the 
vicinity of /ig. The strongest effect is obtained for the 
constant term in these expansions, and this leads to our 
toy example. For a total number n of constant inser- 
tions, we get a contribution to 3S(A;,a;) proportional to 
P^2^hk,uj) - (-l)>^"2^(-fc,-w), where p-^l{k,uj) stands 
for the n-th partial derivative of pp^h(k,ui) with respect 
to LO. Using the fact that p2,i{k,Lo) and pi^2i—k,—Lo) 
behave as (w — /zq)^ for lo close to po, we notice that 
a single self-energy insertion in the sunrise graph pre- 
serves the property that 5S](fc,Li;) vanishes for lo = pQ. 
More generally, this implies that '^'S{k,LO — po) — 
up to third order in perturbation theory. However, a 
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Fermi surface deformation already occurs usually for the 
simplest Hartree and Fock graphs, which are first order 
in the coupling strength. In the standard perturbation 
scheme, the dressed chemical potential and the dressed 
Fermi surface are determined by solving the infinite set 
of equations: 



(C2) 



with the constraint that the total volume of the Fermi 
surface does not change as interactions are switched on. 
We therefore see that, in this scheme, there is no reason 
for which 3S(fcF,w = /i) should vanish. This is rather 
unsatisfactory on physical grounds, since it would imply 
a finite life-time for particle-like excitations lying just on 
the dressed Fermi surface. 



APPENDIX D: FIELD-THEORETICAL RG 

This Appendix is devoted to a detailed derivation of 
the field-theoretical RG equations which we have gath- 
ered in Appendix ^ The two main considerations we 
wish to stress here are: i) the choice of external mo- 
menta in the renormalization prescriptions, which has to 
be adapted to the Fermi surface shape, and ii) the use of 
the logarithmic approximation. 



1. Motivation and general idea 

In the usual "field-theoretical" RG, the high-energy 
Hamiltonian is given and fixed, and one parametrizes the 
theory by low-energy values of proper Green's functions 
(for instance the interaction vertices) at a typical scale v. 
Requiring that all these theories at different energy scales 
should correspond to one and the same high-energy the- 
ory yields RG flows, when v is varied. This approach 
is physically natural, because it is based on the calcula- 
tion of low-energy observables. Furthermore, as we will 
see, it allows for a study of crossovers between high and 
low-energy regimes. 

One of its limitations is that it requires renormal- 
izable interactions (i.e. the existence of a continuous 
limit). But it is not really a severe drawback, since non- 
renormalizable interactions are expected to be irrelevant 
(by power counting) in the low-energy limit, which is the 
most interesting to us. Note that the renormalizability 
constraint disappears in RG schemes based on Wilson's 
idea of gradual mode elimination. Several groups have 
recently implemented Wilson's approach to the RG, ex- 
pressed via the Polchinski equation,-2*2, or its one-particle 
irreducible version."*'"^^ Although these equations are ex- 
act, they are quite complicated, since effective interac- 
tions involving an arbitrary number of particles are gen- 
erated along the RG fiow. Any numerical computation 
therefore requires drastic truncations in the effective ac- 
tion. By contrast, the field-theory approach involves 



only a much smaller set of effective or running couplings, 
which is a good feature for practical implementations. 



2. Renormalization of the interactions 

First of all, we have to define the renormalized cou- 
plings. The two corresponding Green's functions (in real 
space for the direction parallel to the chains) are: 

-(0|T [^k,i+sA^4)Kj-sJX3) (Dl) 

X<,Xp'(^2)<,,.,(Xl)]|0), 

x^[,xp'(^2)<,,,.,(^i)]|0). 

In the above equations, |0) is the interacting ground- 
state, T is the time ordering operator, and X is a short- 
hand notation for {t,x). The renormalized couplings are 
the values of the amputated one-particle irreducible parts 
of these Green's functions, divided by i. In fact, the 
charge and spin couplings are the coefficients obtained 
from the Fourier transform of G^, factor of It-^t'^p,p' 
and (Tt.t' ■ crp,p' respectively. The Umklapp coupling 
Us{I, J) is defined in the same way from , as the co- 
efficient of It r'lp p' (the one in front of Ip t'^t p' being 

-C/j_/_5(/, J)). ' 

The set of external frequencies is chosen to be the same 
for all types of couplings. We have decided to take: 



3zy 

-3 = y, 



(D3) 



and by energy conservation we have of course u}4 = 
LUi + UJ2 — ^^3 — I' is the typical energy scale of 

the interaction process, and is the quantity to be varied 
to get RG flows. 

It is a bit more difficult to choose the external mo- 
menta, because of the warping of the Fermi surface. Our 
choice has been dictated by a few natural requirements. 
First the symmetries of the Fermi surface should be re- 
spected, as the right-left symmetry, the up-down sym- 
metry (i.e. ky ^ ~ky, in terms of the original trans- 
verse momenta). Interactions processes for which it is 
possible to choose all external momenta on the Fermi 
surface, should be computed for this special choice, be- 
cause it would otherwise mean the introduction of a spu- 
rious energy scale. Let us first consider G^*-''-'(/, J). It 
is possible to choose fci — kpj, k2 = —kp^j, k^ — 
—kp^j-s, and k^ — fepj^^, only if momentum con- 
servation fcp,/ — fcp.j = kpj+s — ^F,,7-i5 is respected. 
In general, this will not be possible, for we will have 
Aks{I,J) = {kFj+s + kF,j)- {kFj + kF,j^s) ^0. Notice 
that up to a minus sign and a factor of 2, Aks{I, J) is 
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simply the generalization of Afcp in the two-chain model. 
It is then natural to split this quantity equally among the 
four momenta. One can check that the following choice 
fulfills all the conditions we have mentioned: 



0.1 



fci 

fc2 



Aks{I,J) 



AksiI,J) 



(D4) 



, AksiI,J) \ 
kF,j-s H -. I , and 



— kpj+s — 



Aks{I,J) 



For the Umklapps, the choice of external momenta is 
dictated by the same requirements. The equivalent of 
Aks{I, J) is now Afcf (/, J) = 2n-{kFj + kF.j + kFj+s + 
kF,.j-s), and the natural choice of momenta reads: 



kFj + 
- ( kFT + 



fci 
k2 

kz = kF,J-5 + 
ki = kFj+5 



^k^{I,J) 
4 

Ak^{I,J) 



(D5) 



4 

AkY{I,J) 



and 



The next task is to draw all possible Feynman dia- 
grams, and compute them. One can then establish the 
field theoretical RG flows, requiring the couplings mea- 
sured at two different scales should correspond to the 
same high-energy theory. The couplings' flow equations 
are given in Appendix IE II and are obtained in the one- 
loop approximation. As argued in Sec. ^1 the most 
interesting effects connected to Fermi surface deforma- 
tion appear at the two-loop level for the single-electron 
propagator. For the sake of simplicity, we shall use a 
hybrid scheme, involving a one-loop approximation for 
the couplings and a two-loop approximation for the elec- 
tronic self-energy. In the case where all the couplings 
remain weak, it is reasonable to keep only the domi- 
nant term in the corresponding flow equation. If by 
contrast couplings have a tendency to grow and become 
large at low energies, experience from the Kondo problem 
suggests adding the subleading terms to the couplings' 
flow does not provide a better physical picture. For the 
Kondo problem, the two-loop approximation predicts an 
intermediate coupling fixed pointj^^*^ whereas the low- 
energy physics corresponds to an infinite coupling fixed 
pointf?-'^^ 

We shall not give any technical detail on the derivation 
of these couplings' fiow equations which is standardt^ 
The main new feature is the use of special sets of exter- 
nal momenta, described in Eqs. ljD4|) and (jD5|) above. 
However, it is worth focusing on the / function that ap- 
pears in these equations (the / function is defined by 
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FIG. 38: The dashed curve is the real part of the coupling 
g satisfying I1D7L with initial condition g{Q) = 0.1. Here 
(5/Ao — 10~^, so that the coupling goes to zero singularly 
for t ~ 4.6. The solid line represents the solution of the 
approximate equation where / is replaced by /. The inserted 
plot is a zoom on the end of the flows. 



f{t = ln(Ao/i/), (5) = 1 a u ^ \6\ and otherwise, see 
Appendix IE In fact, the "true" RG equations do not 
involve this function, but rather: 

1 



/(i = ln(AoH,<5) = - 



iri 



177 



(D6) 

It is obvious that these / diverge for j/ = |(5|, so that some 
couplings will diverge or vanish singularly at the scales 
given in ljE4|l to IIE12p . Though this physically signals 
the crossing of the characteristic scales, it is practically 
unpleasant for the numerical simulations. Furthermore, 
if we did not work at zero temperature, the energy scale 
given by the temperature T would suppress these diver- 
gences. Notice that the presence of irj factors also implies 
that the couplings will not remain real. 

For all these reasons, it is thus natural to try to find 
a way to get rid of these singular behaviors. This can 
be achieved by replacing / by a function that extends 
its asymptotic behaviors (for v ^ \S\ and ly ^ \d\) up 
to ly = \S\. This is exactly what the function / does. 
The quality of the approximation can be checked on a 
very simple flow equation, for which one knows the exact 
solution: 



1 



1 



5 + b] 



11] 



g\v). (D7) 



This is simply a RPA-like flow, for only one coupling g. 
We will not study this in detail, but we show the good 
quality of the approximate solution for a positive initial 
couphng, on Fig. |2H1 

The approximation we use is in fact nothing but 
a logarithmic approximation. Indeed, / appears in 
the flow equations, after we have differentiated ln[(i^ — 
5 + 'wi){v + 5 — i77)/(2wFAo)^] factors (with respect to 
the scale v)^ coming from the Feynman graphs' log- 
arithmic divergences. Changing / in / just amounts 
to replacing this logarithm by the approximation 
21n[Max(i/, |5|)/(2wfAo)]. 
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FIG. 39: Sunrise diagram involving two G interactions. 



3. Renormalization of the propagator 



For reasons explained at the end of Sec. IIIIB II the 
one-loop self-energy correction does not have much influ- 
ence on the Fermi surface deformation for the quasi ID 
Hubbard systems considered in this paper. Therefore, 
we will only focus on the two-loop sunrise diagram (the 
Kohn-Luttinger diagram being equal to zero). We only 
show the Feynman diagram here, with all the informa- 
tion about the internal lines, for the "G^" contribution, 
on Fig. 1211 There is of course a "f/^" contribution, and 
the spin algebra has to be taken into account. In order 
to simplify the expressions, we have computed the self- 
energy at fc = fcp,/, and not for a general momentum, 
since we have decided not to take Fermi velocity renor- 
malization into account. We give the expression of the 
self-energy in Aopendix lE 21 From this, the first thing to 
do is to compute the counterterms, so that the inverse 
propagators vanish on the dressed Fermi surface. Focus- 
ing on second order terms only, this amounts to require: 



V/e {1,. 



■,N}, 



fi) - 5fi^^' + vpdkY' = 0, (D8) 



(2) 



together with the conservation of the total particle num- 
ber, is found by summing all these equations: 
(5^(2) = _ ^Rj{k ^ kFj,LU = n)] /N. We let the 

reader check that it is in fact sufficient to take the 
Umklapp contribution in this last equation, because the 
G contribution sums to zero. This is true because of 
the following properties: Aka{I,J) = — Afc_Q,(J, /), 
Gl^'\l,J) = G'}'^{J,I) and G'S'J {I + a, J ~ a) = 

Gl^^\j — a, I + a), so that in the sum over /, J and a, 
the terms (/, J, a) and (J, /, —a) will cancel each other. 

Once the chemical potential is known, the Fermi mo- 
menta counterterms can be found: 

1 



6k\- 



(2) 



VF 



.(2 



(D9) 



In order to save space, we will not give the full expressions 
of these counterterms. 

From all this, we can deduce a dressed propagator, 
as we did in Eq. 14()|) for the two-chain model. As in 
Eq. (|46() . the result is still divergent when Aq is sent to 
infinity. This simply means that the counterterms are 
not sufiicient. Something more is needed, and as is well 
known, this is wave function renormalization. The renor- 
malized (R,I) propagator is defined as usual by: 

G^^f\k,u;) = Z,,jG^j{k,u), (DIO) 

and the wave function renormalization factor Zr,./ is 
found by imposing the following renormalization pre- 
scription: 

-1 



G 



'^j {k ^ kpj, uj ^ fi + u) = 

ZRjG^j{k ^ kpj, uj = ii + y) = v. (Dll) 

We have not written which variables Z^j depends on, 
in order to make the equations lighter, but it should be 
clear that it is a function of the couplings at scale of the 
dressed Fermi momenta, of v and of Aq. The calculation 
of the renormalized propagator is achieved thanks to the 
standard observation that Eq. IjDlip implies: 



(fc,c^;.9;i^,Ao) = (D12) 

{k,uj;g';iy',Ao), with 
Znj{g;iy, Aq) 



ZR,,/(g>',Ao) 
ZR,i{g;iy, Ao) 



ZR.i{g{g;i',i^',AQ);i^',AQ)' ^^^^"^ 

In these equations, g is a shorthand notation for all the 
couplings, at scale v, and g' is the same at scale v' . In the 
last equation, the function g is relating the value of the 
couplings at two different scales by g' = 'g[g;v,v' , Kq). 
Finally, the RG flow equations for the functions 
are found by differentiating the multiplicative relation 
"PKjig'iV^v" .Kq) = LpB^j{g]v,u' ,KQ)LpYij{g']v' ,v" .Aq), 
and one obtains: 

d 

■^mjig; i", , Ao) = <^R,/(5; V, v' , Ao) (D14) 
ov' 

d 

-Q^^'R ji'Sig'^; ^ Ao); v\ v" , Ao) 

Notice that the importance of the functions lay in the 
close link between these and the renormalized propaga- 
tor, that one can deduce from Ea. l|D12|) and the renor- 
malization prescription of Eq. IjDllI) : 



G 



R,L 



{k = kF,uj; g"; v, Aq) = uj <^R.L(g"; w, Ao). 

(D15) 

For the N chains we are interested in, the flow equations 
of p-Rj are given in Appendix EH Eq. ||ET4|| . The flow 
equations of can be checked to be exactly the same 
(this is due to the Left-Right symmetry of the system we 
study). 
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APPENDIX E: RG FLOW EQUATIONS 
1. Flows of the couplings 

The flow equations for the couphngs that are given be- 
low are the field-theoretical RG equations, obtained after 
the general analysis of Appendix IdI has been performed. 
In these, the RG time t is given hy t = ln(Ao/;/) (Aq 
being the ultra-violet cut-off and the typical energy 
scale of the interaction), and the running couplings are 
low-energy interactions G{v). 

The high-energy flow equations of the couplings in the 



cut-off scaling scheme are easily deduced from the ones in 
the field theoretical version. They are in fact the same, 
with the sole difference that the / functions appearing 
in the fiows are to be replaced by 1 (which is the value 
these / functions take at high energies). Notice however 
that the quantities entering the fiow equations now have 
different physical meanings. The cut-off scaling time is 
t = ln(Ao/A) (A being the running cut-off), and the cou- 
plings are running bare couplings Gb(A). 

The field theoretical RG flow equations for the charge, 
spin and Umklapp couplings are given below: 



Q 

|/ (t, 2t>FifP|^5(/, J)) [g^(/, J + a - 5)GUM + a,J) + SGUI, J + a~ 6)0^,^^(1 + a, J) 

-/ (t, 2z;fXPp (/, J)) [g^J/, J)GL J/ + a,J-a)- 3G^J/, J)G^,_J/ + a, J - a) 
+f (t, 2vFK^!i (/, J)) [c/a(/, J + a- 5)Us-M + J) + Uj-isil, J + a- 5)Uj-i-s{I + a, J) 

-^C/5_a(/ + a, J)Uj-i-s{I, J + a- 5)- hjj-i-s{I + a, J)Uc,{I, J + a - S) 



(El) 



a 

I / (t, 2vpK^'^,{I, J)) (/, J + a - <5)G^,_„ (/ + a, J) + G^„ (/, J + a - 5)G^_„ (/ + a, J) 

+GI {I,J + a~ d)Gl_^ iI + a,J) 
+/ (t, 2t;FifPP (/, J)) [2G^(/, J)GL J/ + a,J-a)- Gl{I, J)GL J/ + a,J-a) 
-Gl,{I,J)G%_M + a,J~a) 



(E2) 



+/ (i, 2vFK'^li{I, J)) C/,7-/_5(/, J + a- 5)Uj-i-s[I + a,J)- -Us-a{I + «, J)Uj^i^s{I, J + a - 5) 



'-Uj^i^siI + a,J)U^{I,J + a-S) 
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a 

/ (t, 2«fX^.^i (/, J)) [Gi„ ( J + a - ,5, /) - Gi„ ( J + a ~ ,5, /)] Us-a. [I + a, J) 

+f (i, 2vfK^^^ (/, J)) [G^,_„ (/ + a, J) - G^,_„ (/ + a, J)] C/„ (/, J + a - 5) 
-/ (i, 2t;Fi^Gt'3(^^ J)) [2Gf_ + (5, /)C/j_„_5(«, J)] (E3) 
-/ (i, 2t;Fi^^;y^(/, J)) [2G}_„_5(a, J)C//-„(a + <5, /)] 

+/ (i, 2^;FXGf (/, J)) [CLJI + S + a,I)+ 3G'_JI + 5 + a, I)] Us{I + a, J) 
+f {t, 2vFK^;i\l, J)) [G}_Ja - <5, J) + 3G}_Ja - 6, J)] Us{I, a) I. 



In the above three flow equations, we have used the 
following notations: 

KZ iI^J)-\ i^ks {I, J)- 2 Afc„ (/, J)] (E4) 
(^"^) = i [^^* (^"^) - 2 Afc„ (/, J + a - <5)] (E5) 
^) = J [^ks{I, J) - 2Afc^(/, J + a - <5)] (E6) 

kZ'^' iI^J)^\ [^kf (/, J) - 2 Afc„ (/, J + « - 5)] (E7) 
^a;f ^) = J [AA:y(/, J) - 2Afc„_5(J, / + «)] (E8) 
^) = J [AA:y(/, J) - 2Afc„_,(/, a + 5)] (E9) 



I 

^^^'U, J)^\ [AA;F(/, J) - 2Afc„+5-j(J, a)] (ElO) 
^) = J [AA;y(/, J) - 2Afc„(/, / + a + <5)] (Ell) 
^Sfa. J)^l J) ~ 2Afc„.j(J, a - d)] (E12) 

We refer the reader to Sec. ID 21 for the definitions of 
Aks{I,J) and Afc^(/, J). The / function is defined as 
follows: /(t = ln(Ao/i^), 5) = liiv \5\ and otherwise. 

2. Renormalization of the propagator 

The two-loop self-energy has the following expression: 



J, a 



2 [G^(/, J)G"_„(/ + a, J-a) +3G^(/, J)G'_„(/ + a, J-a)] 



(E13) 



X + vpAkail, J) In 



+ [/„(/, J) [2[/a(/, J) - Uj-I-ail, J)] (i^ + vpAk^il, J)) 



In 



(2i;FAo) 



|z/^-(z;FAfc^(/,J))^ 
(2i;fAo)2 



The functions relating renormalized propagators at in the case of N chains reads: 
two different scales satisfy the general Ea. (|D14|) . which 
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J, a 



2 [Gl{I, J)G"_„(/ + a,J-a) + 3G^(/, J)G''_^[I + a,J~a) 



1 + ^-E^ML^n.] /(.,.pAfc„(I, J)) - ^-I^ML:^iiu,v,AkM, J)) 



+U4I, J)[2UaiI, J) - Uj-I-41, J)] 



(E14) 



where the function I is defined by l{i'^5) — ln(|i^/(5|) if 
h' ^ I (5 1 and = otherwise. This function is, as /, a 
logarithmic approximation of a more complex function, 
diverging at scale \6\. From the definition of (pj (see 
Eg. ID13)l . it is clear that if the initial high-energy quasi- 



particle weight is equal to one, then one simply has : 
Zj{t) = l/^i{t). ^ 

The flow equation for the running Fermi surface in the 
cut-off scaling scheme is also obtained thanks to the self- 
energy, and we find: 



]\f2 A^] a 



2N' ^ I 



+ Afc^^°V> J) Ua{I, J) [2C/a(/, J) ~ J)\ 



(E15) 



- ^ y Afc^ (/, J) [/„ (/, J) [2C/„ (/, J) - (/, J) 
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